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Abstract 

Optical imperfections, misalignments, aberrations, and even dust can significantly limit sensitivity 
in high-contrast imaging systems such as coronagraphs. An upstream deformable mirror (DM) in 
the pupil can be used to correct or compensate for these flaws, either to enhance Strehl ratio or 
suppress residual coronagraphic halo. Measurement of the phase and amplitude of the starlight halo 
at the science camera is essential for determining the DM shape that compensates for any non- 
common-path (NCP) wavefront errors. Using DM displacement ripples to create a series of probe and 
anti-halo speckles in the focal plane has been proposed for space-based coronagraphs and successfully 
demonstrated in the lab. We present the theory and first on-sky demonstration of a technique to 
measure the complex halo using the rapidly-changing residual atmospheric speckles at the 6.5 m MMT 
telescope using the Clio mid-IR camera. The AO system's wavefront sensor (WFS) measurements are 
used to estimate the residual wavefront, allowing us to approximately compute the rapidly-evolving 
phase and amplitude of speckle halo. When combined with relatively-short, synchronized science 
camera images, the complex speckle estimates can be used to interferometrically analyze the images, 
leading to an estimate of the static diffraction halo with NCP effects included. In an operational 
system, this information could be collected continuously and used to iteratively correct quasi-static 
NCP errors or suppress imperfect coronagraphic halos. 

Subject headings: iiistrumentatiomadaptive optics, techniques :interferometric, methods:statistical, in- 
strumentation: miscellaneous , techniques : miscellaneous 



1. INTRODUCTION. 

Extrasolar planets (ESPs) are expected to be 10 6 -10 9 
times fainter than their host stars in the thermal infrared. 
Solar system analogues will have planetary orbita l distri- 
butio ns within one arcsecond of their host star (|Nielsenl 
2011), where they are masked by diffracted and scat- 
tered starlight. Without a coronagraph, most ESPs will 
appear to be several decades fainter than the telescope 
diffraction pattern, characterized by its point spread 
function (PSF). A coronagraph suppresses the diffrac- 
tion structures in the PSF halo, but wavefront aber- 
rations, alignment errors, and dust or other transmis- 
sion flaws within the telescope scatter starlight into the 
suppressed region, overwhelming the faint science signal 
from the ESP. Ground-based telescopes rely on adaptive 
optics (AO) to correct the atmospheric wavefront distor- 
tions, but even a high degree of correction leaves ran- 
dom, rapidly-changing, residual speckles that are several 
decades brighter than the ESPs. Ideally, long exposures 
should average out the speckle noise as 1/y/t to a level 
where the planet becomes detectable against a smooth 
background. Unfortunately, any unsuppressed halo is co- 
herently mod ulated by the AO speckles, increasing th e 
speckle noise (|Bloemhoj l2004t lAime fc Soummerl 12004). 
and AO processing lag can cause the residual speckles 
to be both brighter and persist long er than they other- 
wise mig ht (lFriedlll99tt lAngel 12003). The most serious 
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features in the scattered starlight are quasi-static speck- 
les ( OSS) caused by non- common-path (NC P) aberra- 
tions (jRacine et al.lll999tlHinklev et al.ll2007l) . which are 
not known a priori, change on timesca les from minutes 
to hours, and average slowly — if at all ([M acintos h et alJ 
12001 ISivaramakrishnan fc Qppenheimerl 120061 )! These 
halo features determine the practical detection limit for 
ESPs. 

A number of techniques have been developed 
to address QSS and the effect of slowly-changing 
PSFs over time. These techni ques include Angu- 
lar Differential Imaging (A PI) (jMarois et alJ 12006k 
ISchneider fc Silverstone!l2003D . Simu l taneous Differentia l 
Imaging (SDL) (IBiller fc Closdl2007t IRacine et al.lfl999l 
iMarois et al. 2000), and Local Optimiza tion of Combined 
Images (LOCI) (jLafreniere et al.ll2007D . all with various 
degrees of success. These techniques all post-process the 
science camera data after observation time, and therefore 
can only achieve incoherent gains. AO systems include 
a deformable mirror (DM) that can dynamically mod- 
ify the incoming wavefront, allowing wavefront correc- 
tion to be introduced pre-detection, with the potential 
for coherent improvements in the signal gain. In an AO 
system, after encountering the DM, part of the starlight 
is diverted into a wavefront sensor (WFS), where the 
wavefront shape is measured and processed by a com- 
puter, sending updated compensating shape commands 
to the DM. The resulting flattened starlight wavefront 
may optionally be used to feed a coronagraph or imaged 
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directly onto a science camera. In either case, NCP aber- 
rations that occur downstream from the AO system are 
only evident in the final science images. This can mani- 
fest as distorted and fainter images or coronagraphically- 
suppressed halos that are spoiled by scattered starlight. 
But if instead of correcting the wavefront from the van- 
tage of the WFS, we can try to improve the science image 
(jMalbet et al.lll995l ) by applying bias offsets to the DM, 
we can obtain coherent gains that are significantly better 
than those achieved post-detection. Suppressing the halo 
also suppresses t he interference or "pinned" speckle noise 
(jBloemhoJ 120041 ) . with a corresponding improvement in 
sensitivity. The problem is determining the correct bi- 
ases to apply. 

Uniquely improving the PSF (or the Strehl ratio) or 
suppressing undesired halo in a coronagraph using in- 
formation from the science camera, requires more than 
just intensity measurements: it requires phase as well. A 
full complex halo measurement can be inverse Fourier- 
transformed to estimate the NCP optical flaws, inter- 
preted under the paradigm of Fourier optics. The Strehl 
ratio can be improved by biasing the DM with the op- 
posite of the computed pupil wavefront. Coronagraphic 
halo suppression can be improved by measuring the resid- 
ual halo's complex amplitude and creating a set of "anti- 
speckles" by adding ripp les to the DM to suppress them 
(|Codona fc An gel 2004). It is possible to build a mod- 
ified Lyot coronagraph t hat uses the blocked starlight 
to for m a reference beam ([Angel 2003|; iCodona fc Ange] 
2004) in a focal plane interferometer. However, mea- 
suring the halo with the science camera requires that the 
reference beam be included, spoiling the image unless ex- 
posures are multiplexed between interferometer and sci- 
ence. The two modes will have have similar exposure 
times if we are to measure faint static speckles that are 
of the same level of brightness as a planet, and the cor- 
responding loss of science exposure time leads to a loss 
in sensitivity. 

Another approach is to use the upstream DM to cre- 
ate weak probe speckles in the focal plane with control- 
lable phase. The speckles coherently add to the under- 
lying halo and interfere, modulating the intensity which 
is then recorded by the s cience camera (|Borde fc Traubl 
[200llGive'on et al.ll200l . By using three or more probe 
speckle phases it is possible to measure the complex am- 
plitude of the static halo under the probe speckle. This 
method works very well and is the basis for the extraordi- 
nary laboratory halo suppression result s at the JPL High- 
Cont rast Imaging Testbed (HCIT) (|Trauger fc Traubl 
120071) . designed to simulate a very stable space-based en- 
vironment. Implementing this on the ground has to also 
contend with the presence of residual AO speckles, as 
well as a much less stable environment where at least the 
quasi-static s peckles are changing on the timescale of a 
few minutes (|Hinklev et"all l2007t [Martinez et"aTll2012fl . 
forcing the the probe speckles to be brighter in order to 
speed the measurements. Regardless of how the complex 
halo measurements are made, the same anti-halo servo al- 
gorithms can be applied or, if the the entire complex halo 
is available, Strehl ratio improvements can be achieved 
at the science camera. The downside of using a speckle 
probing method is the same as the coronagraphic focal 
plane interferometer: it puts probe light into the science 
image, forcing science and wavefront sensing to be mul- 



tiplexed, reducing sensitivity. 

In this paper we present and demonstrate a new tech- 
nique for measuring the complex halo in the focal plane 
using only the residual AO speckles. Since the speck- 
les are already present in the field of view, no extra 
light is introduced during the measurement, hence there 
is no need to reduce science exposure time. Since the 
AO speckles vary rapidly, measurements of both quasi- 
static speckles and other halo structure can be obtained 
fast enough to implement focal plane servos with up- 
date rates of several times per minute. Our method 
does not attempt to control the residual AO speckles, 
but only uses the WFS telemetry to monitor and char- 
acterize them as they change. The method requires no 
new hardware and should be implementable on virtually 
any AO system. The only requirements are that the sci- 
ence camera be capable of sufficiently short exposures (< 
30 ms) to adequately "freeze" the speckles, and that the 
camera and WFS data be accurately synchronized. In 
an earlier paper (jCodona et al.l 120081 ). we proposed an 
algorithm called "Phase Sorting Interferometry" (PSI), 
which binned and sorted the pixels based on the com- 
puted speckle phase. That result gave the static halo's 
phase, which is the information needed to properly place 
a DM ripple in an anti-halo servo. In the present paper 
we greatly simplify the analysis into a set of statistical 
formulae that can easily be computed in real time. If 
the halo is estimated over only a limited search area, 
as it would be with a coronagraph, the results can be 
used to suppress the residual halo. If the halo is esti- 
mated more globally, including the PSF core, the result 
can also be used to estimate and correct the NCP aber- 
rations, improving Strehl ratio and PSF quality at the 
science camera. 

Our technique is most closely relat ed to the random 
phase diversity techni ques employed bv lLee et al. I (|1997t ) 
and more recently bv lFrazinl (|2013f ). In those papers the 
focal plane intensity is written as a nonlinear functional 
of the unknown NCP aberrations, as well as the known 
residual AO wavefront error. The NCP error is then es- 
timated from a synchronized set of focal plane images, 
simultaneous WFS measurements, and a nonlinear algo- 
rithm. In contrast, our approach is primarily to create 
the functional equivalent of a focal plane interferometer 
that measures the mean halo field. Once known, the lin- 
ear relationship between the focal plane and pupil plane 
fields can be used to estimate the mean pupil field. If 
desired, the aberrated wavefront may then be estimated 
by computing the complex argument of the mean pupil 
field. Our approach never requires a nonlinear algorithm 
to estimate the pupil wavefront, the phase nonlinearity 
being completely contained within the calculation of the 
pupil field argument. 

In Section 2 we present the theory of interferometry us- 
ing a known-but-random probe field. In section 3 we sim- 
ulate an AO system with NCP aberrations and apply the 
theory to the simulated data set. We demonstrate the ro- 
bustness of the estimated complex halo and derived NCP 
aberrations by adding a wide range of noise in the science 
images. In Section 4 we use Shack-Hartmann WFS data 
from the MMTO 6.5 m telescope and short-exposure im- 
ages from the Clio mid-infrared camera to estimate the 
complex halo. By Fourier transform, we estimate the 
complex pupil field including the effect of subsequent 
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NCP aberrations. Using a 7th magnitude test star, rea- 
sonably consistent NCP aberration estimates were avail- 
able in a few seconds. More stable estimates would re- 
quired longer integrations, but an update to a possible 
NCP-correcting servo would be possible every 30 seconds 
or more. 

2. THEORY 

We begin by considering a set of short-exposure images 
of a star, acquired while the telescope's AO system is de- 
livering diffraction-limited images with a clearly visible 
PSF core. The images consist of a static halo (the diffrac- 
tion pattern of the telescope) and a cloud of rapidly- 
changing speckles. The speckles are significantly fainter 
than they would be be without the AO system, but they 
are still present and contribute scattered starlight out 
to the radius of the seeing disk and beyond. The sci- 
ence camera is unable to directly measure the incident 
electrical field's phase, only the intensity of the static 
halo as it is modulated and coherently interfered with by 
changing speckle cloud. To see the speckles and their in- 
terference clearly, the individual exposures must be short 
enough that the speckles do not change much in phase 
or amplitude. The residual speckles arise from several 
sources: wavefront correction errors from noise in the 
wavefront sensing, limitations or errors in the wavefront 
reconstruction, AO servo errors (e.g. loop gain), fitting 
errors due to the deformable mirror, and processing lag. 
Of these, lag error speckles are the most dynamic since 
they are related to the wind along the telescope's line of 
sight. The speckle cloud typically appears as diffraction- 
limited noise, with individual speckles having widths of 
~ X/D. But the halo also contains speckles and possibly 
deformations and extended regions of scattered starlight 
with much longer timescales. These include the quasi- 
static speckles which confusion-limit sensitivity in high- 
contrast detection. For our discussion, we consider these 
more slowly-changing features to be a part of the "static" 
halo, and by measuring its complex amplitude we can 
provide an adaptive anti-halo servo algorithm with the 
information required to correct aberrations or suppress 
unwanted halo and maximize sensitivity. 

Whatever the cause of the rapidly-changing uncor- 
rected wavefront error, and certainly in the case of lag 
error, the WFS residuals are non-zero and can still be 
measured. While the continuing goal for high-end AO 
systems is to reduce the residual wavefront error to be- 
low the WFS noise floor, we are not there yet. It is 
true that some of the residual error cannot be seen by 
the WFS, but in the current generation of AO systems 
there is always statistically-significant measurable wave- 
front error while the AO system is in closed loop. We will 
assume that we are in that situation. Any unsensed but 
varying residual will appear to us as unexplained noise in 
the image data. Changes that are averaged over during 
the science exposures will appear as a loss of coherence. 
Classic errors such as waffle are not a problem with the 
MMT since the hexapolar actuator pattern is incommen- 
surate with the WFS subapertures. It might be an issue 
when implementing this method on other AO systems 
however. 

Even though the Taylor "frozen flow" hypothesis (i.e. 
the approximation that turbulent motion of the atmo- 
spheric irregularities may be ignored as they are carried 



past the telescope's line-of-sight by the wind) may be 
nearly true for the incident uncorrected wavefront aber- 
rations, the residual AO wavefront error does not behave 
the same way. Consider a typical low-altitude wind that 
carries atmospheric irregularities across the 6.5 m MMT 
aperture in about 0.5 s. High altitude winds are typically 
much faster and cross the aperture in about 100 - 200 ms. 
In either case the irregularities are not likely to signifi- 
cantly rearrange during the pupil crossing. Turbulence 
within the observatory dome or in the vicinity of the 
telescope optics are different and significant rearrange- 
ment is likely. The AO system estimates the wavefront 
error and approximately removes it after a processing 
delay. It does this using a wavefront error-suppressing 
servo with an adjustable gain factor. As a particular 
wavefront pattern crosses the pupil, the AO system it- 
eratively measures the residuals and adds corrections to 
the DM. Even for rapidly-moving wavefront patterns, the 
AO system has many tens to hundreds of iterations to 
suppress it. As a result, the residual pattern is usually 
incrementally altered as it moves across the pupil, be- 
coming progressively more uncorrelated with its earlier 
configuration as it moves. The effect of an incorrect servo 
gain is that any pattern will not be completely removed 
in the next iteration. This repeats even while the ini- 
tial pattern is being carried across the pupil, leaving a 
residual wavefront pattern that exponentially decays in- 
place with a time constant that depends on how far the 
gain is from the "correct" value. Thus the residual wave- 
front error is typically dominated by two or possibly three 
parts: a set of changing wavefront patterns that appear 
and disappear with a time scale dependent on the servo 
gain, and one or two wind-driven patterns that decor- 
relate much faster than turbulent rearrangement would 
suggest. Meanwhile in the focal plane, the wind-driven 
lag error forms a wide plume of speckles that are brighter 
near the star. In addition, as the wavefront errors move 
across the pupil, the speckles' complex phasors rotate 
at a speed which increases linearly with distanc e from 
the s tar projected along the direction of the wind (jAngell 
l2003h . Eventually, the phase wrapping is so fast that the 
speckle phase changes significantly during a single science 
camera exposure, causing the speckle and the static halo 
to appear to lose coherence in that interference between 
them has less contrast. Superimposed on this systematic 
effect is a random evolution of the speckles that depends 
on how long it takes to replace the pupil-plane aberra- 
tions with a statistically unrelated set. This is usually 
a fraction of D/ \\v w i n d\\- These effects, along with the 
general fading of the speckle halo with angular distance 
from the star creates a practical outer radius for using 
speckles as a halo measuring tool. Fortunately, the gen- 
erally brighter and slower speckles nearer the star allow 
us to measure the aberrations that affect telescope or 
coronagraphic performance exactly where improvement 
is needed most. 

In isolation, the phase of a speckle is unobservable with 
the science camera unless the speckle coherently inter- 
feres with the static halo. But the cloud of rapidly- 
changing residual AO speckles do constantly interfere 
with the halo, making an intricate pattern of intensity 
fluctuations that can be captured by the science cam- 
era. Using the WFS telemetry and an idealized model 
for the telescope that ignores NCP errors, we can com- 
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Speckle Phase 

Fig. 1. — Complex halo measurements using MMT data for a 
single pixel in the first Airy ring. The science camera intensity is 
plotted against the speckle phase computed from the WFS mea- 
surements. The speckles vary randomly in phase and amplitude 
while the underlying halo is presumed to be steady. When the 
speckle phase is the same as the halo, the intensity recorded by 
the science camera is increased. When the phase is 180° off from 
the halo, the intensity is reduced. If the speckle amplitude were 
constant, the intensity would vary sinusoidally with speckle phase. 
The randomly-varying amplitude causes the sinusoid to become 
noisy, but the mean behavior is the same. Just judging by eye 
we can see that the underlying halo has a phase of about -2.5 rad 
in this pixel. Statistical analysis allows for much more precise es- 
timation of the halo phase. Intensity scatter vs. intensity swing 
allows us to estimate the underlying halo amplitude. See Sec.|4]for 
details. 

pute the phase and amplitude of the speckles at the sci- 
ence camera. Comparing these calculated speckles with 
the intensity recorded by the camera allows us to use the 
speckles as interferometric probes, enabling estimation 
of the phase and amplitude of the static halo. This re- 
lies on a simple concept: when the speckles are in phase 
with the underlying halo, the intensity is greater, and 
when they are anti-phased the combined intensity is di- 
minished. If we were to plot the observed halo bright- 
ness against the speckle phase (Fig.[T]), the maximum in 
the light curve corresponds to the phase of the halo at 
that point. In addition to the phase, an estimate of the 
static halo's amplitude can be found by comparing vari- 
ous statistics computed from the science images and the 
computed speckles. 

We begin by deriving the interferometry equations 
needed to use an ensemble of known random (in the 
sense of measured but uncontrolled) reference beams to 
estimate a static subject beam. We then apply the re- 
sulting equations to the problem of random speckles and 
a static halo. Since we do not have any a priori knowl- 
edge about the NCP aberrations, our calculations of the 
complex speckle field can only be approximate. How- 
ever, the phase of the speckle peaks and their locations 
are able to be computed rather robustly and the effects 
of the unknown NCP errors will only enter as "speckles 
on the speckles" which will appear as a small amount of 
noise in our measurement. 

2.1. Traditional Interferometry Theory 

Since we have the ability to create speckles in the sci- 
ence camera's focal plane by impressing a ripple on the 



DM, we could simply cycle through three or more ripple 
phases (i.e. by changing the ripple placement relative to 
the pupil) and determine which placement reinforces or 
suppresses the stati c halo under the resultin g speckle (see 
IGive'on et all [20071 ; iKenworthv et al.ll2006[ ). The analy- 
sis is straightforward and equivalent to the common in- 
terferometry technique of using a phase-shifted reference 
beam to probe an unknown subject beam. In that gener- 
alized case the subject beam under test is mixed with a 
phase-shifted reference beam and a camera is used to cap- 
ture the resulting interference pattern or "interferogram." 
Using residual AO speckles is similar to using intention- 
ally created speckles, except that they vary randomly in 
both phase and amplitude and are beyond our control. 
Even though the evolution is random, we can still mea- 
sure the residual wavefront and determine the speckle 
amplitude and phase by calculation. The temporal evo- 
lution of the speckles is only important in that they might 
change during a single science camera exposure, smear- 
ing the parameters that we would like to use in interpret- 
ing the resulting image. Otherwise, the random speckles 
can be thought of as instantaneous samples drawn from 
an ensemble described by statistical moments. Interfer- 
ometry with a random collection of reference beams (or 
speckles) looks different than the conventional analysis 
with carefully selected phases and a constant amplitude, 
but the principles are the same. Therefore, we will start 
by reviewing phase-shift interferometry, but use termi- 
nology that extends cleanly to our random speckle case. 
Both the traditional and generalized theory described be- 
low can be applied to measuring fields in any optical 
plane or physical context, but we will only consider its 
application in the focal plane in this paper. The tradi- 
tional goal is to measure the phase of a constant subject 
beam ^ by mixing it with a constant reference beam if) 
that is stepped through a set of known phase offsets. At 
each phase offset the camera is used to capture an in- 
tensity interferogram that is the magnitude-squared of 
the sum of the subject and phase-shifted reference beam 
fields. In our treatment we consider a single pixel of the 
camera and the fields are complex constants with phases 
4> s = arg'J' and <p r = arg-0 (where for ip = \ijj\ exp{i</>}, 
<fi = arg ip = 3 In ip) . The reference beam is phase-shifted 
by an extra amount 77 £ {?7n}- To solve for the subject 
beam, we need at least three phase shifts in order to 
uniquely determine the phase. Even though everything 
here is constant or preset, we will use angle bracket no- 
tation (•) to indicate averaging over an appropriate or 
indicated set. This will simplify the notational transition 
to randomly-varying reference beams easier. We require 
that the phase shifts are selected such that 

1 N 

( exp {ir, n }) = -J2^ n =0. (1) 

n=l 

If this condition is not met, the non-zero-mean reference 
beam will leave a residual that cannot be distinguished 
from being part of the subject beam. Our analysis will 
therefore produce the correct value biased by the mean 
of the reference beam. We also require that 

(exp{2zr/„}) =0 (2) 

so that we may easily process the interferograms. The in- 
terferograms record the intensity of the interfering fields 



Focal Plane Wavefront Sensing using AO Speckles 



5 



at each of the applied phase shifts, giving 

I n = \* + i>e ir *>\ 2 (3) 

= |*| 2 + \i/)\ 2 + *V*e""'" + l-W" (4) 

= |*| 2 + + 2 IV-I |*| cos(0 s - <f> r - 7] n ). (5) 

The expanded real form in Eq. [5] is the most familiar, 
but we prefer Eq. 0] because the complex algebra is sim- 
pler. We proceed by assigning each interferogram to the 
shifted complex phase of the reference beam and averag- 
ing over the phase shifts 

(J„e"'") = (|*| 2 + |V^| 2 ) <e"'") + + ^Xe 21 "") . 

(6) 

Making use of our constraints Eqs. [T] and [21 we find that 
the mean intensity term averages to zero, and one of the 
two varying interference terms is "stabilized" while the 
other averages to zero by construction. The single sur- 
viving term is the product of the subject and conjugated 
reference beam 

(i n ^) = ^r- (7) 

Taking the imaginary part of the log of Eq. [7J we find 
the subject beam phase to be 

0,=arg(/„e*^> + ^.. (8) 

In many interferometric applications, the only desired 
quantity is the subject beam phase, in which case we are 
done. In cases where the subject beam amplitude is also 
desired, it can be found by simply blocking the reference 
beam and measuring the non-interfering intensity. But 
since we cannot "turn off" the speckles, along with the 
fact that the speckle cloud's shape and amplitude fluctu- 
ates with the wind and other external influences and may 
not be accurately known, we will now find the subject 
beam amplitude assuming the reference beam amplitude 
is unknown. From Eq. [5] we find the average intensity as 

(i) = m 2 +m 2 . (9) 

The magnitude-squared of Eq. [JJ gives us the product of 
beam intensities 

|(l„e^)| 2 = |vI/| 2 |V>| 2 . (10) 

Note that * and ip are interchangeable in both of these 
equations. Using Eqs. [TU] and substituting in Eq. [S] and 
solving for the subject beam intensity, 

|*| 2 = M±i v /(J n ) 2 -4|(/„e^)| 2 , (11) 

we find two solutions. This is due to the subject- 
reference ambiguity built into the statement of the prob- 
lem. Which beam was physically phase shifted in our 
measurements is lost when looking at the intensity, and 
we are left with two solutions: one being the subject 
beam intensity and the other the reference beam inten- 
sity. 

2.2. Interferometry with a known random reference 

beam 

The analysis in Sec. 12.11 is very rigid in its reference 
beam variation, the main purpose being to simplify the 
mathematics. But we can also use a randomly-changing 



reference beam to analyze a static subject field VP in 
each pixel of a set of images in a manner very similar 
to Sec. 12.11 Being "random" in this context means not 
being controlled or predictable. "Random" does not im- 
ply that the reference beam is unknown, since we will 
always assume that at least its phase is known. If we 
know both the reference beam's phase and amplitude 
corresponding to a series of images, the problem is very 
similar to the phase-shifted reference beam case above. 
Unlike that analysis however, which was unable to tell 
the difference between subject and reference beam phase 
and amplitudes, the varying reference beam allows us to 
unambiguously solve for both the subject beam since it 
is constant while the reference beam changes. If all we 
know about the reference beam is its phase, we can still 
estimate the subject beam's phase by using the known 
reference beam phase to stabilize one of the fluctuating 
terms in Eq. 01 Even though the amplitude is changing, 
this allows a reliable statistical estimate of what refer- 
ence phase maximizes the intensity (Fig. [1} . When the 
intensity is maximized, the subject and reference phases 
are the same. If we also know the reference beam am- 
plitudes, or even if we just know the rms speckle am- 
plitude, we can determine the subject beam amplitude. 
Since our reference beam is determined by calculation 
from the WFS measurements and does not depend on, 
for example, the flux of the star (Sec. I4.4[) . we also have 
a scale factor between reference and subject beams that 
needs to be determined from the data's statistics. 

In Sec. 12. 1[ we typically had 3 or 4 images with known 
phase offsets. For the random reference problem, we 
imagine having a synchronized pair of relatively large 
data sets: one of complex reference beam fields and the 
other of the corresponding intensities. We treat each 
synchronized collection of images and fields as a set with 
statistics determined from the available data. We make 
no assumptions about underlying probability distribu- 
tions. 

The total field is the sum of the subject and reference 
fields, with the reference beam changing without the need 
for a separate phase offset 

hiatal = * + V- (12) 

The subject field * is assumed to be constant over the 
data set, while the reference beam tp changes in some 
uncontrolled yet measured or otherwise known way. The 
first assumption is that the reference beam field has a 
zero mean, functionally corresponding to the reference 
phase shift requirement in Eq. Q3 but allowing for varying 
amplitudes, 

W - 0. (13) 

If we were analyzing an ensemble of data, this would 
have some more absolute meaning. But since it is just the 
mean over the available data, perhaps spanning just a few 
seconds, the mean will have some random residual. This 
residual will appear to be part of the static subject beam, 
affecting the derived subject beam as estimation error. 
If the reference beam mean is non-zero, we can subtract 
it from the complex reference values and proceed using 
Eq.O 

*totaI = (tt + <V» + (tf-M). (14) 

The analysis will now yield the biased result (* + (if)}), 
from which the means may be subtracted at the end of 
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the calculation. For simplicity, we will continue to write 
the reference field with the suppressed mean as ip which 
makes Eq.[13]true by construction. In addition to Eq.[T3l 
we note that the mean square (not conjugated) of the 
field becomes negligible in large data sets, 



<^ 2 ) -> o. 



(15) 



This is the analog of Eq. [5] and cannot be independently 
forced like the mean. In our case where the speckles (our 
reference beam) are each the sum of many independent 
patches of approximately-flattened field across the pupil, 
the complex speckle field tends toward a Gaussian dis- 
tribution by the central limit theorem (App. [A}. Note 
that for a zero-mean Gaussian random field, the scatter 
in estimates of (ip 2 )\ oc where N is the num- 

ber of statistically-independent speckles in each data set 
(App. [E]). The same is true of any moment (ip n (ip*) m ) 
where n ^ m. We only require this general result for one 
more term, where the scatter in estimates of (tpip^ip) — > 
as 1/VjV which is likely true so long as the speckle in- 
tensity and phase are uncorrelated (see App. |B|) . 

We start as before with a simple intensity model, leav- 
ing out any questions of partial coherence or intensity- 
dependent noise. We simply treat the static halo as a 
complex constant coherently interfering with the chang- 
ing but known reference field. The resulting intensity is 
given by 



J=|* + V>| 

i= |*| 2 + |^| 2 + *v* + *>■ 



(16) 
(17) 



Taking the average of Eq. [T71 and invoking (ijj) — 0, we 
find the subject beam intensity is 



1*1 



(i) m 



$r — $ 



speckles • 



(18) 



To facilitate our use of these results in later sections, we 
have also written the result in terms of physical measur- 
ables: = (I) is the mean PSF recorded by the science 
camera and ^ spe ckies = (V-'V'*) is the average intensity 
of the halo of speckles if they were recorded on their 
own. Note that while the intensity expressed in Eq. [T7] 
is real, each of the conjugated cross terms *?/>* and ^>*if> 
rotate in opposite directions and individually have zero 
means. We can stabilize the first of the fluctuating terms 
in Eq. [T7| by multiplying through by ip, maintaining a 
stable phase while causing the other term to move twice 
as fast in phase. Upon averaging, the stabilized term 
survives and the others drop out giving 



(19) 



In writing this equation, we also made use of our assump- 
tion that (\ip\ 2 ip) —> even though it will not be precisely 
zero for a shorter data set, but will be in a larger data 

set. Since (IV'I 2 ) € R ; we can immediately determine the 
phase of the subject beam as 



arg * = arg (ipl) , 



(20) 



which is the robust result we described in iCodona et al.l 
(2008). Since the reference beam intensity must be non- 
zero in order to make any measurement at all, we can 



simply write the subject field as 



* 



speckles 



(21) 



2.3. Including an unknown scale factor 



The subject beam field * is the static diffraction halo 
in the science camera focal plane. The reference beam 
field ip is the speckle field caused by the changing resid- 
ual wavefront errors after the AO system has "corrected" 
the atmospherically-distorted incident field. The inten- 
sity / measured by the science camera is proportional 
to the magnitude-squared of the sum of the static halo 
and speckle fields, but with the additional factors of in- 
cident flux and exposure time, etc. The speckle field is 
computed from WFS slopes with possible (modal) gain 
errors in the wavefront reconstruction, and generic as- 
sumptions about amplitude. Converting these fields to 
actual images require all of the aforementioned factors 
such as starlight flux. All of these influences can be in- 
cluded in an unknown scale factor g that may also be 
a function of focal plane position. The relationships be- 
tween *, ip, and / depend on getting this scale factor 
right, so we will add it to the model and estimate it from 
the data. We introduce the scale factor into Eq. [TBI with 
the modification 



I=\V + g^ 



(22) 



where g G K since we assume our phase estimates are 
accurate. Following our derivation as before, we find the 
mean intensity 



*/ = 1*1 + 9 $ speckles 

and the stabilized intensity 

(ipl) = 0* $ speckles- 



(23) 



(24) 

To determine the scale factor before we know *, we ad- 
just it until the intensity variance based on the calculated 
speckles matches that of the actual intensity images. Us- 
ing Eq. [55]to compute a] = ((I — (I)) 2 ), we find 



*1 



2g 2 |*| 2 * 



speckles 



4 2 
y ^speckles 



(25) 



where a 2 peckles = ({\ip\ 2 - (\t/j\ 2 )) 2 ) is the speckle-only 
intensity variance without halo interference. Using the 



magnitude-squared of Eq. 
sured quantities in Eq. [21; 



a 



speckles 



M]to replace g 2 |*| 
5J we find 

2|(^)| 2 

* speckles °~ speckles 



with mea- 



i/4 



(26) 



Using this result in Eq. [24] allows us to write the halo as 



* = 



9^ speckles 



(27) 



This is the diffraction halo at the science camera, includ- 
ing any non-common-path aberrations. While it is obvi- 
ous that the intensity imaged by the science camera must 
first have any background noise removed, it is important 
to note that the background noise variance should also 
be estimated and removed from aj before using Eq. 1271 
We will also find that the region near the center of the 
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halo (very low spatial frequencies) tends to be biased low 
due to physical reasons not included in our isolated pixel 
analysis. Since the result is still useful without increas- 
ing the complexity of the model to eliminate the bias, we 
will use these results as they stand. 

2.4. Visibility 

The traditional interference fringe "visibility" measures 
the intensity modulation caused by changing the ref- 
erence beam phase normalized by the mean intensity 
(|Born & Wolfl 119991 Sec. 7.5). It is usually written in 
terms of the max and min intensities as 



V 



In 



(28) 



Since the standard deviation of the cosine is y/2, we can 
use Eq. [5] to rewrite the conventional definition of visi- 
bility as 

V = (29) 

In our random reference case, the appropriate portion of 
the intensity variance is the term in Eq. [25] that is caused 
by interference (2g 2 \ty\ 2 speckles) rather than intrinsic 
variations in the reference beam intensity (g 4 <J 2 pec ki e s) '■ 
This allows us to write 



V 



2gj\Pj y/Q speckles 
(I) ' 



Using Eq. [23 to replace \ty\, and the notational associa- 
tion of (I) = $/, the formula for the visibility becomes 



V = 



2 |^/>| 



speckles 



(30) 



3. SIMULATION 



In this section we demonstrate our technique with a 
realistic simulation based on the MMTO 6.5 m telescope 
and AO system, with some arbitrary NCP errors added. 
Using the simulation results, we apply our random refer- 
ence interferometric formulae to estimate the static com- 
plex halo and informally explore the effect of measure- 
ment noise. Finally, we Fourier transform the complex 
halo back into the pupil plane to estimate the NCP aber- 
rations. We find that this method underestimates the 
very lowest spatial frequencies of the pupil field, requir- 
ing a simple correction to recover the proper static wave- 
front aberration. The results are adequate for building 
either an anti-halo servo for a coronagraph or improving 
the Strehl ratio at the science camera by compensating 
for NCP aberrations. Finally, we discuss the degrada- 
tion of results caused by speckle halo evolution during 
the science camera exposures and how to estimate its 
effect from real data. 

We used a simple lagged spatial filtering model to simu- 
late the post- AO wavefront at the MMT. We first synthe- 
sized a 2048 x 2048 grid with 4 cm spacing of Kolmogorov 
wavefront displacements scaled to give a Fried length r 
of 15 cm in V band and an outer scale of 30 m. The effect 
of AO correction was simulated by applying a Gaussian 
smoothing to the wavefront and subtracting the result 
from the original wavefront. This approximates the effect 




Fig. 2. — (a) The simulated non-common-path error masked by 
the pupil. The modeled aberration includes trefoil and a Gaussian 
indentation. The peak-to-valley aberration is 725 nm. (b) The 
NCP error added to a realization of the residual AO residual WFE. 



of applying a mode-limited reconstructor to the DM. The 
width of the smoothing kernel was adjusted to make the 
rms residual wavefront fitting error (WFE) 250 nm. The 
effect of wind and processing lag was introduced by later- 
ally shifting the smoothed correction by Sx — v w i n dTi a g- 
For convenience, the lateral shift was selected to be one 
grid pixel of 4 cm. A total of 5 phase screens were synthe- 
sized and a set of residual wavefronts was selected from 
each by using a sub-grid of size D = 6.5 m at 0.85-D spac- 
ings giving a 14 x 14 sample of slightly-overlapping wave- 
fronts. Each scintillation-free pupil field instance was 
computed using exp {2Tti5z/\} and multiplied by a 6.5 m 
pupil mask II(x) with a 10% central obstruction. The 
wavelength A was set to 3.8 /jm (L' band). In addition 
to the residual AO wavefront error, the camera images 
were presumed to include an additional non-common- 
path (NCP) wavefront error described by Zernike trefoil 
(Sz = Z\ (2r/D, 9) A/16) and an additional Gaussian in- 
dentation centered on xo = (—1,2) (meters) with depth 
A/8 (max mirror figure error of 237.5 nm) with a 1/e 
radius of 1 m (Fig. [5]). The PSFs were computed using 
2-D FFTs padded to 2048 x 2048 samples and taking the 
magnitude squared to find the intensity. The resulting 
images were converted to array counts assuming a mean 
total of 157,000 counts per exposure (typical of a magni- 
tude 7.5 star at L' band, bandwidth of 0.65 fim, MMT 
aperture, 50% QE, and a throughput of 67%). Gaussian 
noise was added to the images to simulate the variance 
caused by different amounts of thermal and sky back- 
ground noise. The complex speckle halo was computed 
without the NCP aberrations, but with the same pupil 
mask. This takes the place of the WFS-based halo cal- 
culation in the actual observations. The intentional ne- 
glect of the NCP WFE in the halo calculation emulates 
the fact that we cannot see the downstream NCP er- 
rors when working from the WFS measurements. This 
"bootstrapping" technique of the algorithm should have 
minimal effect on the speckle phases so long as the Strehl 
ratio is high enough to have a clear peak above the speck- 
les. (We tested the effect of this approximation by in- 
cluding the NCP errors in the speckle calculation, with 
an effect on the estimated wavefront error of 8 nm rms 
or about 0.013 radians rms.) The resulting complex halo 
has speckles with the ideal PSF's sidelobes, as opposed to 
the actual speckle sidelobes which would show the trefoil, 
etc. The computed speckles are otherwise in the correct 
locations and have the correct phase relationship to the 
PSF core. After subtracting off the mean complex halo, 
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Fig. 3. — The mean speckle halo intensity averaged over the 
simulated data set. The residual AO fitting error is isotropic, while 
the lag error driven by the wind {v w i nc i T lag = 4 cm) is responsible 
for the horizontal plume of speckles near the center. 

the remaining speckles have mostly the correct complex 
values near their maximum amplitudes, where the effect 
of the simplified sidelobes can be ignored. 

In a calculation of this sort, it is possible to know and 
control the piston component of the wavefront, result- 
ing in absolute phase control or to have a global phase 
reference. This is not the case in reality however. The 
wavefront sensor only measures wavefront slopes and is 
therefore blind to constant phase offsets. The mean in- 
tensity distribution of the speckles, $ speckles , has power 
and phase fluctuations all the way in to the center of the 
PSF. We multiply the complex halo by a unit-amplitude 
phasor to adjust the halo's peak value to zero phase and 
reference the speckles to the current halo peak, allowing 
them to be compared. Ideally, we would want to refer- 
ence the halo phase to the core of the static halo, but it 
is only practical to use the phase of the full computed 
complex halo. This will slightly bias the reference phase 
and introduce a fluctuating phase error to all points in 
the halo. This phase error is zero-mean and becomes less 
significant at higher Strehl ratios like those in our simu- 
lation. Our phase referencing procedure does not change 
either the speckles' intensity or variance, which is why 
we correctly see speckle power all the way to the cen- 
ter of the PSF (Fig. [3]). When we process real data, we 
also remove tip-tilt from the WFS measurements and the 
science images, which has the effect of further removing 
speckle power to the next order about the center of the 
PSFs. Again, this is not detrimental to the measurement 
of away from the center, but it will bias the lowest spa- 
tial frequencies and become an issue when we attempt to 
estimate the pupil field from our focal plane results. 

Once the complex halo was computed for all pupil 
wavefronts and all of the peaks were referenced to zero 
phase, the complex mean was computed and subtracted, 
leaving only the complex speckles -0. The mean speckle 
intensity is shown in Fig. [3J The figure clearly shows 
the broad fitting error cloud and a horizontal plume of 
speckles caused by lag error and the wind. There is no 
evidence of the main diffraction pattern in the speckle 
halo, since the appearance of aberrations in the simula- 
tion is independent of position relative to the pupil. To 
explore the robustness of the result, zero-mean Gaussian 
noise was added to the simulated CCD images to simu- 




Fig. 4. — Simulated science camera images: single frame (left) 
and data set averages (right) . The images are constructed from the 
computed PSFs with photon noise and added background noise. 
The starlight is set to an average of 157,000 counts per exposure, 
with added background noise sigmas of (a) 0, (b) 10, (c) 100, and 
(d) 1000 counts per pixel. The PSF images are 15X/D on each side 
with 0.25A/D pixels. The simulation dataset included 980 images. 
The PSFs on the left are individual 20 ms images while those on 
the right are averaged over the 980 images in the entire data set (a 
total exposure time of 19.6 s). 

late observations of fainter target stars. The introduced 
noise sigmas are 0, 10, 100, and 1000 counts per pixel. 

In order to use the random reference interferometry 
results, we need to compute various statistics from the 
simulated fields and images. The mean speckle field is 
zero by construction, leaving any actual non-zero means 
folded into the static halo as estimation error. The sim- 
ulation generates two data cubes: a set of simulated sci- 
ence camera images with all aberrations included and ac- 
tual counts per pixel, and a data cube of complex halos 
that, in reality, would have been computed from the WFS 
measurements and a simplified model of the telescope 
without the NCP aberrations. The mean of the complex 
halo data cube is an estimate of the static halo, but with- 
out the NCP distortions we wish to measure. This mean 
is subtracted from the complex halo cube leaving only 
the complex speckle field. While the computed speckle 
phases and relative amplitudes values are meaningful, the 
units are arbitrary and cannot be compared directly with 
the photon counts in the science camera image cube. So 
we include a scale factor g (first introduced in Eq. 12121) 
as a constant that depends on the flux and takes care 
of the units. The average of the science images gives $/ 
(Fig. H|), the average of the magnitude-squared speckle 
fields (the "speckle intensities") gives Q speckles (Fig. 
the variance of the science images is a\ (Fig. [5|) and the 
variance of the speckle intensities is u 2 speckles . Since the 
two datasets are computed from the same grids, they 
are automatically commensurate and properly aligned. 
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Fig. 5. — The variance of the simulated science camera images 
<r| displayed on a 4-decade logarithmic scale. This figure only in- 
cludes intensity variations caused by the changing speckles, with no 
background or photon noise. The science image variance includes 
the prod uct of the speckle halo and the intensity of the static halo 
(Eq. I25t . There is a reduction in the variance at the center as a 
consequence of the anti-Hermitian symmetry of speckles and is not 
included in our analysis model. 




Fig. 6. — The computed static halo and visibility at noise levels 
listed in Fig. [4] The leftmost column is the base-10 log of the static 
halo amplitude. The levels are normalized to a point in the first 
Airy ring. The center column is the computed phase of the static 
halo, and the rightmost column is the visibility (lighter is better). 
Note the missing intensity at the very center of the halo in the 
left column images. This is due to the anti-Hermitian symmetry 
of the speckles, an effect that is not included in our interferometry 
analysis. 



We multiply the science image and complex speckle field 
data cubes element- by- element and average over the 960 
frames to compute (tpl) from which we can derive the 
static complex halo if? using Eq. [57] (Fig. \§§ . Note that 
the static halo has a lower amplitude than it should at 
the very center of the function. This is caused by a sub- 
tle difference between our interferometric analysis and 
the symmetries of complex speckles. We assumed that 
an arbitrary random complex speckle could be added to 
the static halo anywhere during the interferometric anal- 
ysis. This is true everywhere in the focal plane except 



at the center of the star image. The reason is that a 
small amplitude wavefront ripple in the pupil plane not 
only diffracts starlight from the PSF core into a speckle 
at the expected position, but creates another speckle on 
the opposite side of the star. Referenced to the starlight 
in the PSF core, the two speckles have different phases 
such that the real part of the field is antisymmetric, 
while the imaginary part is symmetric. This is called 
anti-Hermitian symmetry and is a well-known feature of 
faint speckles. Each speckle is able to take on an arbi- 
trary complex value, subject to energy constraints, and 
the speckle on the other side of the star will follow ac- 
cording to its anti-Hermitian symmetry. This non-local 
correlation is not an issue in our analysis because the 
two correlated speckles are usually so far apart from each 
other that they have no mutual effect. However, a very 
low spatial frequency aberration, or even small residual 
tip-tilts, cause speckle pairs which appear close enough 
together that they can significantly interfere. When they 
do, the antisymmetric real part tends to cancel and the 
sum becomes only due to the symmetric imaginary part. 
Since we have constructed the phase reference to be the 
PSF core, ideally of the static halo center, the center of 
the static halo is real. Adding a small speckle near the 
core approaches purely imaginary and has the effect of 
a phase shift. This does not change the intensity of the 
PSF core, only the phase, therefore the intensity variance 
near the core should drop relative to where the speckle 
pairs don't overlap (Fig. [3]). 

If we use our results in an anti-halo servo designed to 
suppress unwanted residual halo in a coronagraph, the 
low spatial frequency power deficit in Eq. [57] does not 
affect us. However, if we wish to estimate the wavefront 
in the pupil plane, as we might for increasing the Strehl 
ratio in the science images by correcting the NCP aberra- 
tions, or for computing a more accurate numerical PSF 
for post-detection processing of the science images, we 
would like to recover the low spatial frequencies as well 
as the higher ones. We can still do this by reconsider- 
ing how we compute \*f?\ 2 . Instead of using Eq. [24] to 
compute \& , we compute the halo phase arg \& , which is 
what we found in Eq. [20] We then estimate the static 
halo intensity from Eq. [53] which does not suffer from 
the anti-Hermitian speckle effect. To satisfy the require- 
ment that \*f?\ £ K and |\&| 2 > 0, we use a max function, 
giving a better estimate of if? as 



# = i/max{$ ; - g 2 ^ speckles, 0}' 



,'iarg(V>-0 



(31) 



We computed \& using this equation and Fourier trans- 
formed back into the pupil plane to compare with the 
initial NCP aberrations. The results for a set of noise 
levels is shown in Fig. [Jj For each noise level, the pupil 
field amplitude is shown on the left and the phase con- 
verted back to nanometers of wavefront displacement is 
shown on the right, which should be compared with the 
introduced NCP aberrations. The spatial resolution is 
limited because the SNR drops off with increasing radius. 
As the noise level is increased, the estimated pupil field 
and wavefront becomes progressively smoother. The sim- 
ulation is somewhat unrealistic in that the speckles were 
not derived from a Nyquist-limited WFS measurements 
and the wavefront itself has properties that are statis- 
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Fi g, 7. — Fourier transforming the static halo \t computed from 
Eq. 1311 gives an estimate of the static pupil field T for the same 
noise levels listed in Fig. [4] The left column is the amplitude of the 
estimated pupil field and the right column shows the pupil phase 
scaled to wavefront displacement in nm. As the noise level in- 
creases, high angles are lost first, leaving a progressively smoothed 
pupil wavefront estimate. 

tically "cleaner" than real data with residual wavefront 
errors that are not tied to the pupil or actuator locations. 
Nevertheless, the results validate the technique and show 
us what to expect with real data. 

3.1. Effects of optical bandwidth and exposure time 

In the preceding discussions, we assumed that the 
evolving halo field was monochromatic and instanta- 
neously sampled in time, with no blurring or phase av- 
eraging due to speckle evolution during the exposure. 
However, these effects can be important in practice. The 
incident wavefront at a single point in the pupil experi- 
ences an rms phase change of 1 radian after a time To de- 
fined by the phase structure function D ^(v W i n dTo) = 1 , 
where \ w i n d is the characteristic wind velocity projected 
onto the pupil plane. This time scale determines how 
fast WFS exposures must be made as well as actuator 
and AO servo bandwidths. However, the speckles do 
not depend only on the field at a point, but the field 
across the entire pupil. Moving the pattern by a small 
distance, or even possibly by ro, leaves most of the uncor- 
rected incident wavefront unchanged, simply translated. 
This causes speckles to remain in the same place, but 
with a phase shift depending on the location of the given 
speckle in the focal plane relative to the wind. AO resid- 
ual speckles are more complicated because they are the 
result of the initial aberrations and the spatial and tem- 
poral limitations of the AO system. 



For atmospheres with visible ro values of ~ 15 cm and 
winds of ~ 20 m/s, tq > 2 ms. An AO wavefront sensor 
typically has readout speeds of 500 - 2000 frames/s, al- 
lowing the wavefront to be estimated at that frame rate. 
Limited photon flux degrades the wavefront estimate ac- 
curacy with increasing frame rate, as does decreasing star 
brigh tness or increasing background noise (Sandl er et al.l 
I1994D . The MMT AO system is currently at the low end 
of the frame rate range, but still fast enough to capture 
the low spatial frequency wavefront evolution without 
significant averaging. Science cameras are usually used 
with much slower frame rates, although capabilities of 30 
to 50 fps are not uncommon. Broader optical filter band- 
widths pass more light, enabling us to make better use 
of higher frame rates. However, the broader filter bands 
also increase chromatic effects that can lead to confusion 
in the analysis. 

Our MMT L-Band images have optical bandwidths of 
about 20% (0.7 /mi bandwidth centered on 3.7 /mi). This 
introduces chromatic radial smearing where both diffrac- 
tion halo and speckle features appear at points propor- 
tionally farther from the star with increasing wavelength. 
A practical operational constraint is to not allow speckles 
and other halo features to radially blur into each other, 
beyond which they may become confusion-limited. The 
limits implied by this constraint can be estimated by not- 
ing that the typical monochromatic speckle size and sep- 
aration are both ~ Xq/D (Ao is a characteristic wave- 
length in the detected band) , regardless of their location 
relative to the star. At a radius 9, the radial smearing is 
(6X/Xq)6. The non-overlap constraint gives us the con- 
dition 

9 (32) 

The MMT Shack-Hartmann WFS uses a 12 x 12 sub- 
aperture array centered on the pupil, which means the 
maximum measurable unaliased spatial frequency is 6 
cycles/1), beyond which we cannot compute the speckles. 
This corresponds to a fractional bandwidth of about 17%, 
which is reasonably well matched to the L band filter. 
Well inside the radial blurring limit, we are progressively 
more free to ignore the effects of ba ndwidth and sim ply 
use monochromatic Fourier optics (jGoodmanl 119951 ) at 
the characteristic wavelength Ao . 

Depending on the brightness of the star compared to 
the background noise and the shape of the speckle halo, 
the SNR of our measurement varies widely across the fo- 
cal plane. A long science camera exposure will smooth 
the intensity fluctuations caused by the speckles beating 
against the static halo, reducing their contrast and mak- 
ing them harder to detect against the background noise 
and fluctuations. The WFS data, converted to complex 
speckles, also has to be integrated to the science cam- 
era's frame rate in order to compare the complex speckles 
with the images. Our metric for signal is the "visibility" 
(Eq. I3"0"f , but that alone does not express the fluctuating 
signal relative to the background noise. The visibility 
is the phase-dependent intensity variation relative to the 
mean intensity. The proper metric for comparing the 
intensity variations, from which we derive all other in- 
formation, and the background noise is 

SNR^-^ = 2 lf> I . (33) 

Q-noise <&noise \l $ speckles 
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The average intensities, $ no ise and $ speckles, are unaf- 
fected by the exposure time since they are sums of pos- 
itive contributions. The impact of exposure comes only 
from the KV'-OI factor, which we can analyze. Treat- 
ing I(t) and ip(t) to be instantaneous measurements, we 
write the effect of longer camera exposures as an integral 
average of duration T centered on t = 0, and use an en- 
semble average to estimate the effect on | (tp I) | . In Eq. [33] 
we replace / with an integrated / to represent the time 
exposure 

T /2 

I(t)^~[ I(t + t')dt', 
1 J -Til 

and 

m^^f at+t')dt'. 

1 J-T/2 

Multiplying the down-sampled measurements and en- 
semble averaging, the stabilized intensity becomes 

(i>(t)i(t)) -> ^ jj ^ (ip(t + h)i(t + 1 2 )) dhdt 2 . 

Referring back to our simple model for the intensity, 
Eq. 1221 we expand and multiply through by the field, 
performing the ensemble average on each term. Based 
on our assumption that the speckle field is ultimately a 
zero-mean Gaussian random process, the only term that 
survives averaging is I(t) — gty?p*(t). Thus 

mm) fS jf_* 2 r(* a - t 1 )dt 1 dt 2 , (34) 

where T(t) is the complex speckle field's mutual coher- 
ence function (MCF) 

T(t 1 ,t 2 ) = m 1 )iP*(t 2 )). (35) 

The MCF becomes time-invariant, depending only on the 
time difference, T(t2 — t\) with the symmetry T(r) = 
r*(— t), so long as parameters such as wind, ro, and 
AO performance remain constant. We can normalize the 
performance by the assuming that the camera exposures 
match the WFS, which gives the reference value of (tpl) 
as 

(ll>I) -> . 9 *r(0) = speckles ■ (36) 

We divide Eq.[M]by Eq.JSSlto estimate the effect of expo- 
sure time on SNR. Changing time integration variables 
to sums and differences and using the MCF symmetry 
allows us to write the effect of exposure on SNR as 

We will use this result with the actual MMT data in 
Sec. PI 

The residual speckle halo $ S peckies(9) = r(ii,ii;0) 
consists of a broad fitting error halo along with lag er- 
ror speckle plumes in the apparent projected direction 
of the wind. If there are multiple wind streams along 
the line of sight, they will each contribute their own 
speckle plume. Both fitting and lag error phase pat- 
terns are carried across the pupil by their respective 



winds, with characteristic effects expressed in the result- 
ing speckles. For simplicity, we will consider only one 
dominant wind stream with velocity projected onto the 
pupil plane of ~v W i n d- Fitting error has minimal power 
in spatial scales larger than £ao ~ D / V N modes and the 
corresponding halo is therefore dark within X/£ao, or 
at least relatively constant, depending on details of the 
AO system. Processing lag error contributes a resid- 
ual wavefront that is proportional to the gradient of the 
uncorrected wavefront dotted into the wind shift after 
the lag (i.e. Tiag^wind • V</>o) since the wavefront cor- 
rection servo makes the same error in every iteration. 
This leads to a plume of speckles in the direction of 
the wind that becomes brighter as we look closer to the 
star. Depending on the details of the AO system, the fit- 
ting error changes completely by the time the wind has 
carried the wavefront by the AO correction scale £ao- 
Therefore, the fitting error speckles decorrelate after a 
timescale of l A0 j \w winA \ < T fitUng < D/ \\v wincl \\ (a 
fresh breeze of 20 m/s at the MMT would give timescales 
of 43 ms< t fating ^325 ms). Lag error speckles are still 
affected by larger spatial scales and can last much longer 
(App.|Aj). These are certainly completely decorrelated by 
the time the wind has carried the turbulent pattern by 
the outer scale, T wind ^ L outer / \\v wind \\ (possible on the 
order of a second or more at the MMT) . The outer scale 
in all layers of the astronomical AO problem are often 
considered to be less than 30 m, but may be longer in spe- 
cial cases. The decorrelation timescales described here 
vary across the speckle halo around the star, and only 
real data can tell the actual behavior. But we can say a 
few things that ought to be robust statements. The fit- 
ting error speckles should decorrelate fairly rapidly, while 
lag error speckles should last significantly longer. 

The same is not true of the translation effect, which 
is systematic and highly dependent on position. Trans- 
lation of the residual wavefront by an amount <5x causes 
any speckles at angular position 6 = n/k, k = 2ir/\, to 
undergo a phase shift of 8(f) = —k9 ■ <5x. (This ignores the 
overall change in coherence due to slightly different areas 
of the wavefront being visible through the pupil at differ- 
ent times, which is contained in the previously discussed 
portions.) Therefore, a wind will cause a linear phase 
shift Scf> — kO ■ v W i n dSt. This affects all of the speckles 
together systematically, and shows up in the MCF as 

r(r) =r Q (T)e 2tkTe - v — (38) 

where r = t 2 — t\ . 

We can now estimate the effect of science camera expo- 
sure time. We will consider only a single wind stream of 
~10 m/s, noting that for a moderate frame rate of ~ 30 
fps, a jet stream contribution with a speed of ~60 m/s will 
be highly averaged over virtually all of the focal plane, 
appearing as incoherent noise. For fitting error, the slow 
change in the MCF due to different configurations being 
present within the pupil has a timescale of ~650 ms but 
likely much shorter due to the AO iterations and other 
effects. Lag error speckles within the control radius but 
outside of the PSF core arise from spatial scales smaller 
than the pupil diameter D and are most likely uncorre- 
cted much beyond the pupil. Therefore they too have 
a maximum lifetime of hundreds of milliseconds, which 
is common. The systematic phase rotation caused by 
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the wind increases most rapidly in the direction of the 
speckle plume and is not seriously detrimental until the 
phase wrapping from the time center of the exposure to 
the endpoints is 7r/2, beyond which negative contribu- 
tions are included. This means that the phase rotation 
timescale, providing us with a conservative exposure time 
limit, is when it > k 10 • v wind \t sci or 



< 



D 



\/D 



(39) 



Using the WFS Nyquist radius for the MMT 12 x 12 
WFS along either axis, ||0|| = 6X/D, and a wind of 10 
m/s, we find that our exposure limit is about 54 ms. As- 
suming continuous exposures, that corresponds to about 
18.5 fps. Referring back to the effect of decorrelation on 
the SNR, Eq. 133 we can see that our limit only drops 
the SNR by 64% at the limiting radius. Using a higher 
frame rate camera is therefore not extremely important 
for sensitivity, but would be for an increased measure- 
ment radius or for higher wind speeds. 

4. OBSERVATIONS AND ANALYSIS 

To measure the complex halo at the science camera, 
we require three capabilities: (1) acquire short exposures 
with the science camera, (2) AO wavefront sensor data, 
and (3) provide a mechanism for synchronizing the re- 
sulting data sets together. Our mid-IR camera is capable 
of reading out small regions of the sensor at frame rates 
in excess of 30 Hz, and our high-speed AO WFS and 
subsequent processing software has an engineering diag- 
nostic mode capable of saving the full system telemetry, 
including raw WFS pixels and computed slopes. Tight 
synchronization between the science camera and the AO 
system is not provided by the MMT and is added using 
system handshaking and logging modifications. We de- 
scribe our solutions to these generic problems here, as 
well as other difficulties to do with the engineering state 
of the MMT at the time of our observations. 

4.1. Facilities 
4.1.1. The MMT A O System 

The MMT Adaptive Optics system (jWildi et al.ll2002D 
is the world's first telescope to use a deformable sec- 
ondary mirror to provide wavefront correction. This ap- 
proach minimizes the number of warm optical surfaces 
between the sky and the science camera, greatly reduc- 
ing the thermal emissivity of the telescope, and makes the 
6.5 m MMT aperture competitive wit h larger telescopes 
for t hermal infrared observations ({Lloyd-Har t et al.l 
2000). The 640 mm diameter deformable secondary con- 
sists of a thin shell mirror 2 mm thick, supported above 
a Zerodur reference body by a fixed central hub, and de- 
formed by 336 actuators in a modified hexapolar pattern 
(Fig. [5]). The actuators provide non-contact forces to 
the shell via electromagnetic voice coils acting on mag- 
nets attached to its inner surface. The gap between the 
shell and the reference body is measured by capacitive 
sensors at each actuator, and is actively maintained by a 
40 kHz servo in a dedicated mirror controller. The typ- 
ical time for the shell to reach a desired position is less 
than 1 ms. 




Fig. 8. — The MMT's deformable secondary mirror with an 
oblique view of the 336 voice coil actuators and their modified 
hexapolar placement. 



The MMT AO processin g is performed by a real- 
time Linux computer system ( Vaithceswa ran et al.ll2008l) 
which reads the WFS camera, computes the required 
wavefront updates, and sends updated commands to 
the deformable secondary controller, all synchronously 
clocked by the WFS frame rate. The 12x12 s ubaper- 
ture Shack-Hartmann WFS (jMcguire et al.l H999) is nor- 
mally operated at 527 Hz. The AO computer calcu- 
lates the wavefront slopes, reconstructs an estimate of 
the wavefront, applies a conservative modal filter, and 
sends the updated actuator commands to the deformable 
secondary controller. 

The MMT's Shack-Hartmann WFS uses a 12 x 12 
lenslet array to image the starlight within each subaper- 
ture onto the center of a binned 2x2 pixel "quad-cell." 
Local wavefront slopes cause the image to shift, chang- 
ing the relative amount of starlight entering the vari- 
ous quad-cell pixels. The pixels are exposed, read out, 
summed, differenced, and mapped through a lookup ta- 
ble matched to the seeing level, yielding estimates of the 
x and y wavefront slopes over each of the 144 subaper- 
tures (|Hardvlll998D . The pixel sums and differences are 
normalized by the sum of the quad-cell counts, plus a 
small bias term. The bias ensures that unilluminated 
quadcells generate zero slope estimates. The resulting 
288 x and y slopes are serialized into a single column 
vector (the "slopes vector" 0). The AO system estimates 
the residual wavefront by multiplying the slopes vector 
with a wavefront reconstructor matrix, Wao, resulting 
in an estimate of the residual wavefront error at the ac- 
tuator positions: 5z = Wao®- An AO reconstructor 
typically has a combined legacy of optical measurements 
and analytical processing, characterized ultimately by a 
number of singula r value decomposition (SVD) modes 
(|Brusa et al.l I2003D . The MMT AO reconstructor uses 
lowest-energy mechanical modes of the thin shell as a 
basis set. This restriction means that the wavefront spa- 
tial frequencies are not uniformly corrected within the 
modal control radius of the mirror, implemented as a 
safety precautio n. Subsequent deve lopments with the 
LBT AO system (|Esposito et al.ll2010D and Magellan ob- 
servatories are less constrained. The wavefront correc- 
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Fig. 9. — Locations of the MMT deformable secondary mirror 
actuators (circles). At the time of our measurements, there were 13 
disabled actuators (light gray circles) that float to a point between 
the adjacent actuators. 34 actuators had problematic or failed 
capacitive sensors (dots), but were still able to be driven in open 
loop by the mirror controller. For the 56-mode correction applied 
by the AO control servo, having 10% failed actuators has only a 
marginal effect. 

tion is post-processed using a "modal filter" to redun- 
dantly ensure that damaging stresses are not applied to 
the shell. Once calculated, the residual wavefront esti- 
mates are multiplied by an overall (scalar) gain factor 
and added to the correction already applied to the mir- 
ror. The result is then sent to the mirror controller where 
it is applied with its own high-speed servo control loop. 

There were several issues with the MMT AO system at 
the time of observation. The gains on individual mirror 
actuators had not been recalibrated for several years and 
in some cases had drifted by as much as tens of percent. 
Thirteen out of the 336 actuators were non-functional 
(Fig. [5]). The loss of these actuators did not significantly 
impact low-order correction modes as the DM shell above 
a deactivated actuator floats to a position interpolating 
that of its neighbors. Another 34 actuators had inoper- 
ative capacitive sensors, requiring their positions to be 
set in the high-speed mirror controller by dead reckon- 
ing rather than using the closed-loop control servo. In 
addition to some mount vibrations, our observations in- 
clude significant pointing swings deliberately introduced 
as a signal for calibration of a pyramid wavefront sensor. 
These vibrations have amplitudes in excess of 100 mas 
(i.e. on the order of a science camera image diffraction 
width). 

4.1.2. The Clio mid-infrared Science Camera 

Our science camera was the mid-infrared Clio system 
((Freed et alJ 12004 . iSivanandam et all 12006( 1. Light from 
the deformable secondary directly enters Clio's dewar 
through a tilted dichroic window, forming an //15 image 
at the first focal plane. A re-imaging lens forms another 
pupil plane, usually for bandpass filters and Lyot stops, 
but also where Clio's Apodizin g Phase Plate (APP) is 
located (|Kenworthy et alll2007( l. A final lens images the 
pupil plane onto the focal plane imaging sensor. Clio's 
sensor is a HAWAII-1 HgCdTe array with 18.5 /im square 
pixels cooled to 75.6 K. The detector gain is 4.9 e~/dn 
with a bias of 3700 dn and a saturation level of 55,000 dn, 
giving a full-well capacity of 51,000 dn (250,000 e"). The 



read noise for a single frame is 19 dn (93 e~). The dark 
current is 50 dn/s (245 e~/s) at 75.9 K. Clio's plate scale 
was measured as 0.0299"/pixel or 4.0 pixels/ (X/D) in L 
band. Clio is controlled by a computer running Linux, 
taking exposures and saving FITS image cubes asyn- 
chronously from the AO system. Data is taken using 
a special 54x108 "sub-stamp mode" to achieve a higher 
frame rate. 

4.2. Synchronizing WFS and Science Data 

Since the speckle phases change rapidly, both randomly 
and in a deterministic way depending on position within 
the field and the projected speed and direction of the 
wind, small synchronization errors can cause both ran- 
dom and systematic biases in the measured halo phase. 
Without taking particular care, we might expect the 
phase error to be determined by how much the speckle 
phase changes during the science camera's exposure time. 
However, this can be reduced by carefully synchroniz- 
ing the WFS frames to the start and end of the science 
frames, thereby achieving phase accuracies more in line 
with the the speckle phase change during a single WFS 
frame. This improvement (a factor of 6.5 in this dataset) 
only helps reduce biases. The loss of coherence caused 
by the science camera exposure time is not recoverable 
and still reduces the measurement SNR. However, by re- 
moving systematic effects, averaging will still yield more 
accurate measurements at the higher frame rate. 

All AO system updates occur at the WFS frame rate. 
In an engineering diagnostic mode, the AO host com- 
puter buffers 10,000 frames of engineering data into a set 
of RAM-based circular buffers. Once these buffers are 
filled, a separate processing thread writes them out to 
the AO computer's hard disk. The MMT WFS frames 
are normally taken at a steady rate, but an engineering 
issue resulted in every other WFS frame being dropped 
from the data telemetry, resulting in the system operat- 
ing at half its normal speed. The result was that the nor- 
mally 527 frames/s rate was reduced to approximately 
263.5 frames/s with a slopes file saved every 38 s. To 
avoid timing problems and other possible buffer overrun 
issues, we did not include any science data sets that ran 
across WFS slopes files. These issues do not impact our 
technique. Clio buffers 100 contiguous exposures in RAM 
before saving them to disk. When the RAM-based im- 
age buffer is filled, the camera acquisition halts while the 
images are saved. The non-standard small images were 
acquired at approximately 40.5 frames/s, filling the 100- 
image buffer in about 2.5 s and saved to disk every 3 s. 
Within each Clio data cube, the images are continuous 
sequential exposures, with each exposure ending as the 
next began. 

Both the AO and Clio control computers run Net- 
work Time Protocol (ntp) clients. The resulting filesys- 
tem and internal header timestamps are directly com- 
pared between computers, allowing us to uniquely asso- 
ciate log entries and data files. Using ntp alone does 
not guarantee synchronization to the WFS frame level. 
The AO engineering diagnostics files contain only data 
and no internal timestamps or other useful FITS header 
information. The primary identifying information was 
a unique file name which included a numerical code de- 
rived from internal counters and the time of the first 
frame, recorded to 1 second accuracy. The operating 
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system also encoded timestamps in the hlesystem inodcs 
when the file was written. Since these are volatile upon 
copying and archiving, we saved them to a file after the 
run using the Linux command "Is -It --full-time > 
timestamps . log" which preserved the file creation times 
in ISO format to the accuracy configured in the kernel 
(only 1-second accuracy for both systems). The Clio file 
creation timestamps are also saved as a fallback pro- 
cedure, but the primary record there was a timestamp 
saved in the FITS file header. 

WFS Slopes Measurements (10.000 vectors) 























L 1 1 II 

HI II 























Clio Science Images (100 images) 

time 



Fig. 10. — WFS and science data cubes have asynchronous place- 
ment and non-guaranteed timing. Although clock synchronization, 
network handshaking, timestamps, and logging procedures were 
implemented, a data-based cross correlation procedure was also 
required to ensure proper synchronization. 

In addition to the timestamps, we implemented a sim- 
ple network handshaking protocol to record the WFS 
frame information to higher accuracy. At the beginning 
of a set of science images, The science camera control 
computer sends a UDP packet to the AO host computer. 
Once received by the AO system, this information is writ- 
ten to the AO log file along with the base name of the 
WFS engineering files currently being recorded and the 
current WFS frame number. This information is later 
extracted from the log and used to align the data sets 
to within a few WFS frames. The handshake procedure 
contains an unknown delay from the time Clio transmits 
the UDP packet and the AO system wrote the current 
WFS frame offset to the log file. This correlation gave 
us the mean number of WFS frames per science expo- 
sure, (tgd/twfs), and the lag between the logged frame 
number and the correct value. 

Using the images from a single 100-frame Clio file, we 
determined the location of either the star's peak pixel 
or its centroid. We then interpolated this time series 
based on an assumed value for (t SC i/t w f s } to the WFS 
frame rate, allowing a cross correlation with the mean 
9 X and 9 y time series determined from the appropriate 
subset of the WFS slope file, extended before and af- 
ter by a generous set of additional frames. The concur- 
rent engineering tests were introducing a lot of tip-tilt 
noise, giving us a strong signal for the timing calibra- 
tion. Our exposures were short enough that this did 
not cause us any problems. Even so, the normal tip-tilt 
noise in the AO system would have been adequate for 
this determination. The science camera and the WFS 
are not necessarily aligned, but were in this case. Even 
so, we correlated all four combinations between the two 
inputs to ensure that we had the coordinate mapping 
correct. By varying the presumed ratio of frame rates, 
we determined (t sci /t wfs ) sa 6.48 ± 0.03. Thus, in 100 
continuous science exposures, assuming the center expo- 
sure was correctly placed, the start and end exposures 




Fig. 11. — Clio star images. These 16 sequential images are 
shown with a 3 decade log grayscale. 



have a placement uncertainty of ±1.5 WFS frames or 5.7 
ms. The cross correlation analysis also determined the 
lag from the handshake protocol to be 3 WFS frames. 
Once calibrated, the handshake protocol alone was suffi- 
cient to synchronize the Clio images with the WFS data 
to within a WFS frame on the average. Without the ran- 
dom WFS readout timing problems encountered during 
our test run, and a measured frame rate for the science 
camera, all the exposures would be placed relative to the 
WFS with equal precision, and overall accuracy would be 
determined by the accuracy of the calibrated handshake. 
This is expected to be within a single WFS frame, which 
is the limit of accuracy for the system. As we shall see 
below, the speckle coherence time varies across the field, 
but the minimum values are comparable to the timing 
accuracy errors at the bounds of the exposures. The re- 
sult is that there were some coherence losses in SNR due 
to temporal misalignment between the datasets, but they 
were not serious enough to limit our results. 

4.3. Science Camera PSFs 

The observed star is a 7.4 mag G5 star at zenithal an- 
gles ranging from 13 to 14 degrees (Fig. ITTjl . Clio was 
config ured for direct imaging with a Barr MKO L' band 
filter (fTokuna ga et al J 12002). centered at 3.8 fim and a 
bandpass of 0.7 /im. At the center wavelength, the Clio 
pixels subtend A/4. 0D. The 54 x 54 pixel region read- 
out from the sensor was calibrated with dark and flat 
images, with noisy or bad pixels flagged and replaced 
by the median of a 3 x 3 window centered on the pixel. 
Image vibrations were removed by shifting the PSF cen- 
ter with bilinear interpolation to a fixed location. The 
PSF peak was located using correlation with a Gaussian 
reference peak. The pointing vibration during our obser- 
vations had a swing of approximately 1.5 - 2 A/ D which 
contributed only a small amount of image motion blur 
during any given Clio exposure. This vibration did not 
significantly impact the data reduction process, but did 
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have the effect of dithering over any bad pixels, reducing 
their influence. It was not possible to precisely know at 
observation time where we were in the WFS buffer, so 
no attempt was made to avoid the times bridging slopes 
files. To avoid any possible timing complexities, any Clio 
image cubes that bridged WFS file boundaries were sim- 
ply discarded. 

4.4. Complex Speckles from WFS slopes 

We first use the WFS measurements to estimate the 
varying wavefront displacement in the pupil plane z(x, t), 
where x is the position in the pupil plane. Using the 
science camera filter band's mean wavelength A, we 
compute the phase shift caused by the displacement, 
<p(x, t) = fcz(x, t), where k = lixjX. Ignoring scintil- 
lation, we write a unit amplitude complex field in the 
pupil plane as exp {i<p(x, t)}, pass it t hrough the pupi l 
stop II(x), and then use Fourier optics (|Goodmanlll995f ) 
to find the complex halo in the image plane 

t) = J e 4K X e lv(x '* ) n(x) d 2 x (40) 

where k is spatial frequency in the pupil plane and 6 = 
n/k is the corresponding angular coordinate (in radians) 
measured from the star. Over a given observation period, 
we can make a distinction between the "static halo" and 
"speckles" as 

*(M) = + V'(M) • (41) 

static halo " speckles" 

(A more detailed discussion of the separation of the 
halo into its static and speckled components is given in 
App. El) Basic Fourier optics is a very simple model for 
an imaging system, with no place to properly introduce 
aberrations that occur after the starlight passes the pupil 
without significantly increasing the complexity of the op- 
tical wave model. Aberrations that occur anywhere other 
than the pupil plane (or images thereof), will affect the 
halo and speckles differently depending on where the star 
appears in the camera's field of view. Since the field of 
interest in high-contrast imaging is very small, we only 
introduce negligible errors by treating all aberrations as 
having taken place in the entrance pupil plane, including 
any downstream aberrations not seen by the AO system 
(i.e. NCP aberrations). It is because of this that we may 
measure the complex halo in the focal plane and then use 
Fourier optics to find an equivalent distorted field in the 
pupil plane. Since we do not know the NCP aberrations a 
priori, we must compute the speckle halo without them. 
Depending on the magnitude and nature of the NCP 
errors, this can have a large effect, distorting both the 
static and speckle halos. At higher Strehl ratios (where 
the peak of the diffraction-limited PSF is much brighter 
than the speckles), each individual complex speckle can 
be thought of as a translated, scaled, and phase-shifted 
copy of the speckle-free halo. Therefore, including the 
NCP aberrations would mostly just alter the "speckles on 
the speckles", adding noise to the speckle halo, but hav- 
ing only a minor effect on the brighter speckle peaks. The 
NCP errors also alter the static halo, often in ways which 
cannot be ignored. However, by estimating and removing 
the idealized static halo from our calculation, we are left 



with just the speckles. This justifies why we can ignore 
the NCP aberrations when computing the speckles, al- 
lowing us to bootstrap our interferometric measurement 
of the true static halo. If we need more accuracy, we 
can iterate the process by using the first estimate of the 
NCP aberrations in a second calculation of the complex 
speckle halo. 

As we collect WFS data, we repeatedly carry out the 
above steps to build a "data cube" of complex halos that 
do not include the effects of NCP errors, but do include 
residual AO speckles along with a simplified static halo. 
This static halo is presumed common to all frames of the 
data cube ($ nmffl = ^(k8 nm ,t w ), where 9 nm is the grid 
of pixel coordinates) and can be estimated by averaging 
over the time index. Since we force the time average 
of the speckles to be zero over each data cube, we are 
actually pushing any estimation error onto the derived 
static halo. This is a consequence of processing with 
smaller data cubes, and the resulting individual static 
halo estimates from multiple data cubes, when averaged, 
should lead to a more accurate answer. 

As mentioned earlier, the WFS slopes are related to 
wavefront displacements by a "reconstructor matrix." 
This is essentially an integrator, working on the set of x 
and y slopes returned from the WFS sub apertures. The 
x and y slopes from the 12x12 WFS are serialized into a 
single 288 x 1 "slopes vector", 0, which is integrated into 
a wavefront by multiplying by a "reconstructor matrix" 
W: z = W®. Since the MMT AO wavefront recon- 
structor, W AOi only corrects a limited number of modes 
(currently 56), it is not sufficient for our complex halo 
calculation. If we consider each pair of modes to be the 
equivalent of controlling the real and imaginary parts of 
the complex amplitude of a speckle, each of which has 
a width of roughly X/D, then 56 modes will take us out 
to roughly 3.7X/D from the star, while the WFS is ca- 
pable of at least 6X/D. To allow us to properly recover 
the residual wavefront to the accuracy of the WFS rather 
than being limited by the AO system, we require a bet- 
ter reconstructor, W r residual- Since our reconstructor is 
not used to drive the DM, nor is it intended to be it- 
erated by being placed inside of an AO servo loop, it is 
not as important to be conservative about issues such as 
modal gain, making it easier to derive an adequate recon- 
structor. Also, since the computed displacements are not 
going to be used to drive the actual DM, we can choose 
to reconstruct the wavefront at more conveniently-placed 
locations across the pupil: e.g. on a square grid, instead 
of the actual hexapolar locations of the physical DM. 
But here, to facilitate direct comparison with the AO 
system's data as well as other live reconstructor tests, 
we used the physical actuator locations. For each WFS 
slopes vector, the estimated residual wavefront displace- 
ments at the actuators are given by z = W res iduai® ■ 

We estimated and removed the mean x and y slopes 
from the slopes vector before multiplying by ~W residual-, 
resulting in a vector of wavefront displacements with 
tip-tilt removed. For computing the complex pupil 
field and halo, we interpolated the wavefront displace- 
ment estimates at the actuators to a 4 cm square grid 
using Delauney triangularization and cubic interpola- 
tion, yielding a wavefront displacement z nm at x„ m = 
(n,m)Sx + Xoo- The complex pupil field was computed 
by T nm = n nm exp(27riz„ m /A scje „ ce ), where IT nm is the 
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pupil transmission mask interpolated to the grid coor- 
dinates. The resulting complex mesh was zero-padded 
to N x N and Fourier transformed, giving focal plane 
samples spaced by 5k = 2ir/N5x and angular spacing of 
59 = X SC i ence /N5x. We adjusted N to match the L-band 
Clio plate scale of 0.25X/D. The focal plane field was 
computed using 2-D FFTs and the results kept in a com- 



plex data cube 



$(K x ,K y ,t w f s ). This processing 



was performed for each set of 100 science camera images 
(~ 2.5 s), with the corresponding set of WFS frames se- 
lected and used to compute a halo data cube spanning 
the science camera exposures at the WFS frame rate. 
Since the synchronization between the WFS data and the 
science camera is only accurate to about one frame from 
the center of a 100-exposure image set, but the duration 
of the individual exposures is 6.5 frames, the starting 
frame index is assumed to be on a frame boundary, with 
subsequent exposures placed on the timeline as they fell. 
The full-speed complex halo was then down-sampled to 
the science camera frame rate by summing the complex 
frames, including linearly- weighted end frames according 
to the computed endpoints. Finally, since the actual pis- 
ton and tip-tilt may have changed during a single science 
camera exposure, we performed the complex halo sums 
before normalizing to the peak phase. Since tip-tilt was 
removed in the slopes before computing the halo, motion 
blur was not fully represented, but the image wander is 
much smaller than X/D during any single science camera 
exposure and the error introduced is small. The complex 
halo estimate at the science frame rate was 



where = 
=fr arg{(e 



^)} = o, 



i«6t sci W 



(42) 



Since (</)) 



=£> arg {(e'*'') j = U, core phase referencing was per- 
formed on the complex WFS PSF halo. 

Since the individual complex halo cubes were so short, 
we did not force the mean speckle field to be zero for 
each individual cube, giving a changing estimate for the 
simplified static field used in the calculation. Instead, 
we used an estimate of the ensemble average of the halo, 
(i.e. J/ pl / 2 *E r o as described in App. [Al and used it to es- 
timate the speckle halo over each cube. This is better 
than using the per-cube average since the estimation er- 
ror appears as a constant across the pupil rather than 
randomly textured (Sec. I4.6j) . Our resulting speckle field 
was 

VW = *nm„ " ^ 1/2 *o- (43) 



Note that while this is no longer zero mean on the scale of 
individual image cubes, it should approach zero mean in 
the ergodic limit of many data cubes. The speckle field is 
now able to be directly compared with the science camera 
image cubes as they are synchronized and have the same 
sampling in space and time. 

Before continuing on to compute the static halo, it is 
useful to look at the behavior of the complex speckles 
and their spatial and temporal statistics to understand 
how our result is affected by exposure time, wind, etc. 
Fig. shows an angle-time cut through the center of 
the speckle field in the projected direction of the wind to 
illustrate the lifespan of the speckles, with both ampli- 
tude and phase variations. The speckle amplitude varia- 



50 100 150 200 250 300 350 400 450 

Clio Frame Numbers 

Fig. 12. — A cut through the computed speckle field for five 
non-contiguous Clio image sets. The upper panel shows the speckle 
amplitudes as a function of time while the lower panel is the speckle 
phase relative to the light in the PSF core. The 100 Clio frames 
of each set are contiguous, but not across boundaries at multiples 
of 100. The speckle amplitudes are arbitrarily normalized to the 
brightest value. 



Computed from WFS Actual Clio Image 



Fig. 13. — Star images computed from the residual WFS mea- 
surements (left column) alongside the corresponding Clio images 
(right). The computed images mainly differ from the actual im- 
ages in optical bandwidth and non-common-path (NCP) aberra- 
tions. The computed images are monochromatic and assume that 
the estimated wavefront seen by the WFS is the only source of 
aberration. 

tion timescale is reasonably consistent across the speckle 
cloud, while the speckle phase has a similar pattern su- 
perimposed on a steady phase rate that increases as we 
move away from the star. These speckles coherently add 
to the as-yet unknown static halo, the intensity being 
recorded by the science camera. We can make a visual 
comparison of the science images frame-by-frame to the 

WFS-based images using ^ nmfi , which is based on our 

simplified no-NCP error model (Fig. [T^| . This figure il- 
lustrates if we are indeed computing speckles when and 
where they were actually observed, and the extent of 
any deviations. For each row in Fig.[l3l the PSF derived 
from the WFS slopes is shown on the left, while the cor- 
responding Clio image is on the right. The two images 
exhibit similar behavior, but with the single-wavelength 
computed PSF having more well-defined speckles than 
the actual science images due to the bandwidth of the 
science images. The speckles in the computed images 
may also be more or less prominent than the real im- 
ages since the reconstructor gain may be incorrect as a 
function of position around the star (i.e. a possibly in- 
correct modal gain). Also, since NCP aberrations and 
low-order effects such as defocus are not included in the 
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computed PSFs, there might be noticeable differences be- 
tween computed and real images, as well as complicated 
differences with the speckles. In our present case the 
comparison is reasonably good, providing a sanity check 
that the synchronization and WFS-to-speckle calculation 
was performed correctly. 

The various statistics are straightforward to compute 
from the image and complex speckle data cubes. The 
intensity variance is 
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the mean speckle halo intensity is 
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the variance in the speckle halo's intensity 
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and the stabilized intensity is 



1 N 

w — 1 



(44) 



(45) 



(46) 



(47) 



The scale factor g is computed using Eq.J^Uand the static 
halo "J is given by Eq. [27J 

The wind speed gusted by more than a factor of 2 and 
changed direction by more than 45 degrees during the 
observation, causing clearly-visible effects in the speckle 
field statistics. When the wind was slower (say 5-10 
m/s), the speckle cloud was more isotropic with a form 
characteristic of fitting error, while stronger wind gusts 
enhanced the plume of lag-error speckles. Fig. [T4l shows 
two selected times with slower and faster wind. For 
each case the figure shows the WFS slopes and three de- 
rived speckle statistics. The 12 x 12 sub-aperture Shack- 
Hartmann WFS slopes were organized into 288-element 
vectors and concatenated into a matrix of slopes (up- 
per images). The speckle halo was computed for each 
slopes vector and the three statistical images were aver- 
aged from the 650 complex speckle frames. The lower 
left images are the average speckle power centered on 
the star, showing the morphological change between the 
more isotropic fitting error halo and the lag error speckle 
plume. The middle image is the average amount of phase 
rotation seen during a Clio exposure, computed from the 
two-time MCF (Eq. [35]): <ty(0,t aci ) = arg{r(0,t sci )}. 
The phase rotation is mostly a systematic translation ef- 
fect, while the loss of coherence due to random effects 
is better described by the drop in the magnitude of the 
MCF |r(0,r)/r(0,O)|, shown in the lower right images 
in Fig. [T3] for a lag of 6t w f s . The linear grayscale runs 
from completely incoherent (black) to completely coher- 
ent (white). The results show that while fitting error 
speckles rapidly lose coherence, lag error speckles remain 
coherent longer in the direction of the wind. This is due 
to the well-known phenomenon of the AO system repeat- 
edly making the same lag-induced correction error as the 
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Fig. 14. — The effect of wind speed on the distribution and 
evolution of the speckles. Two 2.5 s WFS data sets are shown, 
with the 2.5 s of WFS slopes vectors organized as a 288 X 650 
matrix. The slopes were reconstructed into residual wavefronts 
and used to compute the complex speckle halo. In the top set of 
images the wind is slower, while in the lower set the wind is gusting 
to more than twice the upper speed. The speckle halo's mean 
intensity is shown on the left. The slower wind case is dominated 
by the relatively isotropic fitting error, while the faster wind shows 
a prominent lag error speckle plume in the (shifted) direction of 
the wind. The average speckle intensity is shown using a square 
root scale. The systematic phase smearing during a Clio exposure 
is is shown in the center image. Finally, the speckle decorrelation 
on the right shows the loss of speckle coherence at a lag roughly 
equal a Clio exposure. 



aberrations are carried across the pupil. 

A more quantitative view of the speckle evolution and 
decorrelation is shown in Fig. [15] As described by 
Eq. [351 P ar t of the MCF describes the speckles' decor- 
relation, while a phasor factor systematically rotates the 
phase depending on the wind's projection onto the se- 
lected point's position relative to the star. The figure 
shows the MCFs, normalized by the zero-lag value (note 
r(0; 0) = $ speckles (kO)) , f° r a single selected point in 
the second Airy ring for all of our speckle data cubes. 
This complex function shows the speckle decorrelation 
as a decrease in the absolute value of T(r) with increas- 
ing time, as well as the systematic phase rotation caused 
by the wind. But also, the figure shows that when the 
wind is faster, the MCF falls in magnitude much more 
slowly. This is due to the AO system repeating the same 
error as it "chases" the wind-driven aberrations across 
the pupil. While this effect increases speckle noise due 
to fewer statistically-independent speckles in a given ex- 
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Fig. 15. — Normalized mutual coherence function r(r)/r(0) for 
a single point in the second Airy ring. Two example cases are 
highlighted, one where the angle between the wind vector and 
the pixel vector is small, (brown with diamonds) and one where 
the angle is significantly larger (blue with circles). For each of 
the WFS-derived speckle data cubes, T(r) was computed at lags 
t : {0, 1, 2, ■ ■ ■ , 20, 25, 30, ■ ■ ■ , 100}t m/s . The speckles always decor- 
relate with time, but wind causes a systematic phase shift in a di- 
rection related to the projection of the wind on the selected point 
in the focal plane. The interesting cases where the angle swerves 
from one sign to the other are caused by more than one wind stream 
competing. The values of the MCF at the end of the Clio exposure 
are shown as filled red circles, showing that for this pixel the phase 
smearing was never a serious issue. The speckle coherence time, 
ro, is defined as when r(r)/r(0) crosses the red dashed circle at a 
radius of e~ x / 2 . The more wind-blown lag error speckles remain 
coherent longer since the AO system repeats the same error as the 
wind carries the aberrations across the pupil. This function and its 
two main beh avio rs (straight decorrelation and systematic phase 
shift, see Eq. I38j l determine th e los s of SNR for various science 
camera exposure times (see Sec. I3.lt . 

posure, it helps us in that longer-lived speckles can be 
imaged with more modest science camera frame rates, 
like 30 fps. Both speckle decorrelation and phase rota- 
tion affect our interferometric measurement's SNR, as we 
will see below. 

In the pre- AO pupil plane, the coherence time is usu- 
ally defined in terms of advection of the wavefront by 
pure Taylor flow. This implies that the natural def- 
inition of coherence time is -D^v^^To) — 1, or as 
is often qu oted, rp = 0.314r /v u ,i rlt j for a Kolmogorov 
wavefront (jGreenwood fc Fried! |l976l ). However, this is 
more properly defined in terms of the temporal MCF, 
T(r) = (ip(x + v wind T)ip* (x)) = exp{-£>0 (v wmci T)/2}, 
where the canonical structure function definition be- 
comes |r(r o )/r(0)| = e- 1 / 2 . This is the definition 
we used to compute the focal plane coherence time in 
Fig. fTCTa). Note that the coherence time is only about 
35-50 ms except in the speckle plume where the coher- 
ence time is much longer. This is the expected behavior 
with a non-predictive AO servo algorithm that uses only 
the most recent WFS measurements to update the DM. 
It is interesting to note in Fig. [THT a) that the lag-error 




Fig. 16. — Average speckle halo characteristic timescales across 
the focal plane (dashed circles at IX/ D spacings). The wind-driven 
speckle plume varied about a direction roughly 60 degrees clockwise 
from the positive y-axis. Image (a) is the characteristic coherence 
time computed using |r(r)/r(0)| = e" 1/2 . The color scale is la- 
beled in milliseconds and runs from 35 ms to 90 ms. The fitting 
error speckles typically decorrelated after 35-50 ms while the lag 
error speckles in the apparent direction of the wind lasted longer. 
Note that the lag-error speckle lifetimes were longer when they 
appeared over an Airy ring, which was not expected. Image (b) 
shows the average time required for the speckle phase to change 
by ±7r/2, extrapolated from the phase swing after 4 WFS frames. 
The color scale runs from 50 to 250 ms. 



speckle coherence time is longer over the Airy rings. This 
was not expected, but may be due to the wind bringing 
in unseen aberrations at the edge of the pupil. More 
study is required. 

The other timescale of concern is how long it takes 
for the systematic phase shift to change by tt/2, after 
which at least some of the interferometric reference beam 
contributions will start to subtract. Note that outside of 
the speckle plume, the phase shift timescale is longer due 
to the oblique translation geometry or irrelevant due to 
the lack of persistent speckles. For speckles that are more 
in-line with the wind flow, the phase change becomes 
more consistent and the lag error speckle plume makes 
the effect dominant. For our data set however, the most 
important effect over the majority of the focal plane is the 
changing residual aberration pattern, not the translation 
by the wind. 

The most useful and practical metric is the additional 
loss of SNR due to longer science camera exposures. As 
expressed in Eq. [33 the SNR is reduced by a weighted 
integral over the MCF at various lags. This can be com- 
puted from our WFS data and is shown in Fig. [T7] The 
figure shows that although the Clio exposure times were 
long enough to show a great deal of coherence loss end- 
to-end, the SNR dropped to only about 60% to 90% of 
an arbitrarily high frame rate camera. This can be re- 
covered everywhere, if required, by integrating about 3 
times longer. 

In each of these figures, there are occasional resonant 
speckles caused by processing lag and wind, along with 
the particular characteristics of the reconstructor, lead 
to enhanced speckles at certain angles. This was an un- 
usual situation, possibly due to the engineering tests be- 
ing performed. However, such phenomena can occur in a 
live AO system and may cause difficulties in the measure- 
ment. In this case, the speckle phase distribution about 
the mean was still sufficient to give a good statistical 
measurement of the static halo. But they are certainly a 
sign that greater care should be taken. 

4.5. Estimating the static halo 
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Fig. 17. — The pixel-by-pixel worst-case loss of SNR due to sci- 
ence camera exposure time, compared to an ideal short-exposure 
SNR. The Clio exposures were 6.5 times longer than the WFS ex- 
posures, so there was some loss of SNR depending on the wind and 
focal plane position. The SNR was only down to about 90% nearer 
to the star, dropping to about 60% in the processing-lag speckle 
plume as the phase wrapping rate increased farther from the star. 
For this bright star and weather circumstances, we would have no 
difficulty analyzing low-order aberrations with even significantly 
longer exposures. 




Fig. 18. — Map of the scale factor g(0) in a neighborhood around 
the star. The scale is base-10 log and the circles are 1—6 X/D. The 
white area on the upper left is due to a dithered flaw in the sensor. 



As derived in Sec. 12. 2[ we can estimate the static halo 
from the Clio and complex speckle data cubes using 
Eq. 1271 In discrete form this is 



w Ynmw ± nmw 



In 



Qnrn 



(48) 



The scale factor is a function of position and is com- 
puted using Eq. [26l shown in Fig. [TH As expected for a 
reconstructor that is not carefully calibrated, the speck- 
les do not have a uniform scale over the working field. 
Including this correction, we now have an estimate of 
the static complex halo as shown in Fig. [T5] Armed 
with this information, we can now use it in an anti-halo 
servo to suppress the h alo over some region of interest 
(jCodona fc Ang el 2004), or we can use the entire halo to 
estimate the wavefront in the pupil plane including the 
NCP error to improve the Strehl ratio and imaging qual- 
ity at the science camera. This use of the measured static 
halo requires more work but is conceptually as simple as 
taking the Fourier transform of the static halo. 
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Fig. 19. — Average complex halo derived from 67 2.5 s Clio image 
cubes (6700 images) and 43,550 WFS measurements. This calcu- 
lation used the presumed ensemble average method for computing 
the speckles and therefore may be missing an additive fraction of 
the ideal complex halo. 



4.6. Estimating the pupil field 

The static halo is th e m ean of the complex halo 
(App. [A] specifically Eq. IAl[) , and therefore it is the in- 
verse Fourier transform of the mean complex pupil field 
T(x, t) under the paradigm of Fourier optics (Eq. IA22j) . 
Our static halo field includes the effects of any non- 
common-path aberrations after the AO system and the 
science camera, which may have been introduced at any 
point in the intervening light path. The inverse Fourier 
transform of the static halo does not give the literal pupil 
field, but the apparent pupil field including the effects of 
the downstream aberrations back-projected to the pupil. 
Since the pupil is the canonical location of the DM where 
we can affect wavefront corrections, we can use this esti- 
mate to determine the wavefront offsets required to mod- 
ify the resulting halo in any desired way. Expressing this 
in terms of the forward Fourier optics equation, Eq. 140) 
we write 



T (x)n(x) 



C *(K)d 2 K 



(49) 



where we have made use of the assumption that the 
residual wavefront phase from the AO correction is a 
zero-mean Gaussian random variable and (exp(z0)) = 
exp(— (0 2 ) /2) = .y 1 / 2 is the square root of the Strehl 
ratio. This gives us an estimate of the pupil field, blurred 
due to the limited spatial frequencies recovered by the in- 
terferometry calculation. The NCP pupil field is ^(x) = 
a(x) exp (ikSz) where the amplitude a(x) G R may in- 
clude transmission effects, as well as back-projections of 
downstream vignetting, etc. The wavefront to be cor- 
rected is given by 



Sz 



2^ 



arg 



(50) 



This simple procedure should be all that is required 
to compute the wavefront correction, but there are some 
caveats and corrections. The pupil field estimate is miss- 
ing some low spatial frequency components, mostly a 
constant. The reasons for this have already been dis- 
cussed in Sec. G2 and are related to our computed phase 
referencing to the PSF core, and the fact out tip-tilt sta- 
bilization reducing the apparent effect of anti-Hermitian 
speckle pairs near the optical axis. We also did not ac- 
count for interference between the very low spatial fre- 
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Fig. 20. — Examples of 2.5 s estimates of the pupil wavefront. 
The scale factor g(0) was computed from the full set of 67 100- 
exposure image cubes and used in the calculation of the static halo 
for each cube. The loss of low spatial frequency power due to the 
inability to sense piston as well as tip-tilt stabilization was approx- 
imately compensated for by adding a constant to the computed 
pupil field. The resulting phase was converted to mirror surface 
heights using z = <£(x)A/47r. The scale is labeled in nanometers. 



quency speckle pairs, causing an error that becomes sig- 
nificant within X/2D of the PSF core. We can calibrate 
this error by applying a known bias to the DM and seeing 
what the measured value is. Since a constant bias would 
affect the measured phase aberration according to 



5z = arg <{ tp bias + 
Zn 



C *(K)d 2 K 



(51) 



we can then adjust the constant ipbias to give a reasonable 
estimate of the actual aberration. 

Since we did not have a calibration for the data pre- 
sented in this paper, we estimated a constant pupil field 
of 50% of the mean computed value. Adding this to the 
computed value for each 2.5 s Clio data cube, we compute 
the NCP aberrations. A selected subset of the wavefront 
errors are shown in Fig. 1201 Since the estimates are gen- 
erated from only 2.5 s duration science camera cubes, 
the recovered NCPs show a significant amount of vari- 
ance, but they also show consistent morphological pat- 
terns. The average of these is shown in Fig. [2TJ The 
peak-to-valley surface error (i.e. half the wavefront er- 
ror) is 1266 nm and the rms error is 206 nm. This result 
is compatible with the science images, and all that re- 
mains is to apply the corrections to the DM and see the 
improvement in the image quality. 

5. DISCUSSION 

We have developed a new interferometric technique for 
focal plane wavefront sensing using the residual starlight 
speckles left behind by an AO system. Our method re- 
quires no extra hardware or optics beyond that used 
in a typical AO system. The rapidly-changing resid- 
ual AO speckles are present in the science camera image 
plane and coherently interfere with the starlight diffrac- 
tion halo. The speckles result from uncontrolled errors 
in the correction of the atmospheric aberrations, result- 
ing in small residual wavefront fluctuations in a pupil 
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Fig. 21. — The MMT's average static pupil field including back- 
projected non-common-path aberrations. The pupil field was bi- 
ased by a constant of 1/2 the rms amplitude to approximately 
compensate for underestimated low-spatial-frequency power. Im- 
age (a) is the intensity, and the 0.1-D central obstruction is clearly 
seen. Image (b) is the average WFE computed as a bias to be 
applied to the DM. The scale is labeled in nanometers. The peak- 
to-valley surface error is 1266 nm and the rms error is 206 nm (an 
rms phase error of A/9.2 at L band). 



plane downstream from the AO system's deformable mir- 
ror (DM). While uncontrolled and unsuppressed by the 
AO system, these residual errors are still monitored by 
the wavefront sensor and are available for other purposes. 
Improved AO technology and algorithms will continu- 
ally reduce the errors and subsequent speckles, but they 
will remain a ubiquitous feature that can be exploited 
for other purposes. Our technique uses the wavefront 
sensor measurements to compute a numerical analogue 
of the complex halo (with both amplitude and phase), 
which is also measured in intensity by the science camera. 
The computed halo is used as a key to interpret short- 
exposure science camera images, allowing the steady por- 
tion of the star's complex halo. Since the intensity mea- 
surements are downstream from the AO system, the our 
complex halo measurement includes information about 
both common and non-common-path aberrations. This 
information is useful both for correcting the wavefront 
before imaging and post-detection image processing. 

We demonstrated the technique using data taken with 
the MMT AO wavefront sensor (WFS) and the Clio sci- 
ence camera in the L band. We synchronized and down- 
sampled the WFS measurements to match the Clio frame 
rate of about 40 frames per second. Using our mathemat- 
ical theory, we were able to compute an estimate of the 
complex halo every 2.5 s from 100 Clio images and 650 
WFS measurements. The resulting complex halo esti- 
mate represents the average field in the focal plane and 
is related to the pupil field by Fourier transform. We 
analyzed the temporal statistics of the speckle halo and 
used them to estimate the effect of slower science camera 
exposures. We found that the impact of exposure time, 
while depending on the details of wind and focal plane 
position, was not serious — especially near the star where 
halo suppression will be most important. 

Although we were not able at this time to demon- 
strate the use of the complex halo to improve the sci- 
ence images, the path forward and potential benefits of 
continuous end-to-end optical quality measurements and 
updating are clear. The Fourier transform of the halo 
provides an estimate of the pupil field, including the 
back-propagated effects of NCP aberrations. By apply- 
ing appropriate DM offsets to compensate for the static 
pupil field's phase, the Strehl ratio seen by the science 
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camera can be improved, providing correspondingly im- 
proved detection sensitivity. Measurements of the com- 
plex halo are even more directly useful for suppressing 
residual diffracted starlight in a flawed or poorly-aligned 
coronagraph by using an antihalo servo, the theory of 
which was not discussed in this paper. 

In future work we plan on developing a system to esti- 
mate the static halo in real time during observations. 
The estimates will then be fed back to the AO sys- 
tem either as WFS or DM biases. For normal (non- 
coronagraphic) imaging, the goal will be to optimize 
Strehl ratio at the science camera. When using a coro- 
nagraph, our goal will be to suppress unwanted halo in 
the search region while simultaneously maximizing Strehl 
ratio. These enhancements will not require the addition 
of any new optics, but will require the addition of an 
outer control loop that affects the operation of the AO 



system. We also plan to use the complex halo measure- 
ments to develop two new image processing algorithms. 
By using the estimated pupil field, we can compute better 
numerical PSFs, which include the effects of both com- 
mon and non-common-path aberrations as well as op- 
tical filter bandwidth. Subtracting the computed PSFs 
from the actual star images will reduce speckle noise. We 
also plan to use the computed PSFs as a key for when 
to keep pixels in the science data cube, including them 
only when the computed halo falls below some statistical 
threshold. This "lucky pixel" type algorithm should be 
capable of increasing sensitivity by a factor of about 2. 
Including both algorithms, the decrease in speckle noise 
can be significant, depending on the raw performance of 
the AO system. For the MMT, we would be hoping to 
achieve an additional sensitivity of 2-4 magnitudes over 
the corrected coronagraph alone. 



APPENDIX 

A. THE STATIC HALO AND SPECKLE CLOUD 

We treat residual AO speckles to be additive random field fluctuations coherently interfering with a static diffraction 
pattern. In this appendix we explore the halo, and consider assumptions and statistics. We show why the high Strehl 
ratio focal plane field is a somewhat attenuated version of the telescope's diffraction pattern, including non-common- 
path aberrations, and a coherent cloud of speckles which are reasonably described as zero-mean Gaussian complex 
field fluctuations. We will consider two cases: the classic seeing-limited situation with no AO wavefront correction as 
a touchstone, and the diffraction-limited case with a static PSF surrounded by a cloud of speckles. 

We start with an incident field T(x) passing through a possibly complex pupil mask n(x), which also includes any 
non-common-path (NCP) aberrations. We will not include scintillation here, but could easily do so. For convenience 
and wavelength independence, we use the spatial frequency k = (k x , n y ) to describe the complex halo field 

*( K ,t) = J e JK X T(x)n(x)d 2 a; (Al) 

where 

T(x) = e^ (x ' t) T (x) 

is the pupil field and T (x) is the pupil field without any of the fluctuating phase aberrations <p(pt,t). Note that our 
halo field is the angular spectrum of plane waves comprising the pupil field, where each k plane wave travels from the 
direction 6 = n/k with k = The point spread function (PSF) is the focal plane intensity distribution resulting 

from an on-axis star, where To(x) = 1. Although non-common-path aberrations occur downstream from the AO 
system and the nominal pupil plane, our small field of interest allows us to treat them as being included in the initial 
pupil field To(x). In our notation, the PSF is the angular power spectrum of plane waves 

*(M) = I*(M)| 2 - (A2) 

The fluctuating phase aberrations are presumed to vary randomly in space and time, with ergodicity allowing us to use 
ensemble and temporal averages interchangeably. While this may be a reasonable assumption for the atmospherically- 
distorted incident field, it is wise to keep an open mind regarding the post-AO statistics. In general, the average PSF 
is 

($(«)) = J d 2 Xl J d 2 x 2 e lK (xi - X2 '> ^e 4(v(xi '* ) - v(x ^ t) ^T (x 1 )TS(x 2 )n(x 1 )n*(x 2 ). (A3) 
The average phase exponential is traditionally simplified using 

<e 4 «) = e-<" 2 > /2 (A4) 

which is valid when q is a zero- mean Gaussian random variable (GRV). If we assume q = <p(x.i) — <p(x 2 ) is a GRV, 
then 

e i(v(xi)- ¥ j(x 2 ))^ _ e -D v (xi,x 2 )/2 ^g-j 

where 



D y (x 1 ,x 2 ) = ((^(x 1 )-^(x 2 )) 2 ) (A6) 
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is the structure function of the pupil field phase. Again, while Eq. IA5I can reasonably be applied to the seeing-limited 
case, the post- AO case may contain patterns and correlations that break this assumption. This may be particularly 
important when using this result with multi-segment telescopes where the segments are being phased with a slower 
servo loop, etc. We cannot address that case here, just provide a warning. The statistics of the uncorrected atmospheric 
phase are also usually considered to be translation independent, while the post- AO statistics may well not be; post- AO 
correction residuals vary depending on position relative to actuators, WFS sub-apertures, pupil edges, etc. A simple 
example is the mean-square phase residual, ^<ys 2 (x)^ = <r 2 , which is generally assumed to be position independent, 
although it probably is not. We note but ignore this for now, assuming that at least the second-order phase statistics 
are translation independent. If so, then 

D y (xi,x 2 ) -> £> v (xi - x 2 ). 

Because the post-AO halo contains lag error speckles, we cannot gene rally assume isotropy since the wind defines a 
preferred direction. Translation independence allows us to simplify Eq. IA31 Make the coordinate change 



a = (xi +x 2 )/2 
/3 = x 2 — Xi 

and carry out the a integration. The result gives the average PSF in terms of a spatially-filtered version of the optical 
transfer function (OTF) 



($(«)) = J e lK - /3 e- D ^/ 2 0(f3)d 2 /3 (A7) 
where the OTF (jWilliams fc Becklundl[2002[ ) is given by 

0(/3) = y"T (a + /3/2)n(a + /3/2)TS(a + /3/2)n*(a + /3/2)d 2 a. (A8) 



Since D v (/3) starts off at zero when (3 — and rises to a maximum value as ||/3|| increases, Eq. IA7I shows how the 
mean PSF loses higher angular frequencies as the rms pupil wavefront phase increases. This expression will be our 
guide in dissecting the halo, both for the intensity PSF and realizations of the complex halo itself. 

The uncorrected structure function has certain common features that we often refer to. As two phase measurement 
points become more distantly separated, the phase values eventually become completely uncorrelated and 

D v (h) -> 2a%. (A9) 

In the post- AO case, this asymptotic limit applies as s — ||s|| Iao where Iao is the correlation scale for the AO 
system, defined below. In the no-AO seeing-limited case, this does not apply until the spacing exceeds at least the 
outer scale Lq, and possibly not until significantly beyond that, depending on definitions. It is not uncommon to 
consider a case, albeit non-physical, where the outer scale and the phase variance are infinite. For simplicity in this 
discussion we will assume that there is an outer scale and (^(x)) = 0. 

The no- AO case is often modeled, or at least described, using the isotropic Kolmogorov power law structure function, 
parametrized by the Fried length ro, 

5/3 



D^s) = 6.88 . (A10) 

This f unct ion actually flattens out past Lq, but for D/ro 3> 1 it will not change this discussion. The blurring factor 
in Eq. IA71 exp {— D v ((3)/2}, has a radius defined by D v ((3)/2 = 1 or j3 = 0.477ro- For our estimates, we will use 
the diameter of the blurring filter as ~ ro. The radius of the OTF is the pupil diameter D , although it may be 
more complicated for a large segmented telescope pupil. If -D/ r o 3> 1, the blurring filter limits the integral in Eq. IA7I 
and we can approximate the OTF with its value at the origin, which is simply the area of the pupil weighted by its 
transmission mask. Therefore in the seeing-limited no- AO case, the average PSF is approximately given by 

($(«)) « 0(0) f e iK - e- D ^/ 2 d 2 /3. (All) 



When the telescope is equipped with a good AO system, the residual phase structure function has a maximum that 
is still relatively small, causing the blurring factor to look like the constant exp{— cr 2 } with a small rise up to the 
value of 1 as f3 —> 0. This no longer limits the integral in Eq. IA7l and reveals a diffraction- limited PSF plus a scattered 
light halo. This effect of this can best be seen by explicitly pulling out the constant term and letting the small rise 
near the origin stand alone 

e -2>„(/3)/2 _ y + _ s^)H{fl) (A12) 

where 

S> = e-^ 2 ) (A13) 
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is the Strehl ratio as usually expressed in Marechal's approximation {Hardy 1998S), and the factor 

-D v (P)/2 _ cp 

H((3) = - — jf (A14) 

goes from a maximum value of -ff = 1 at /3 = to .ff — > for ||/3|| > £ao- Just as the structure function exponential 
limited the integral (Eq. IA7|) in the seeing-limited case, H{0) limits the integral in the scattered halo term when 
D/£ao ^ 1. This gives a high-Strehl approximation for the PSF as 

($(«)) pa y<5>v{n) + O(0)#(k) (A15) 

where 

H(k) = J e lK l3 H(j3)d 2 (3 (A16) 

and 

*o(«) = / e lK/3 0(/3)d 2 /3 (A17) 



is the diffraction-limited PSF (including NCP effects since it includes To(x)). Note that Marechal's approximation is 
always true, but only for the diffraction-limited pa rt of t he PSF. 

So in both the low-Strehl seeing-limi ted case Eq. lAf f l and the high-Strehl post- AO case Eq. lAf 51 we have a spatial 
filter that limits the /3 integral in Eq. IA7I In one case the spatial filter di amet er is r , while in the other it is £aoj 
the width of H((3). Although they were derived somewhat differently, Eq. IA15I actually applies to both the AO and 
no-A O cas es, corresponding to the high-Strehl and low-Strehl cases respectively. 

Eq. IA15I describes the mean high Strehl ratio PSF as the sum of a static PSF plus an averaged speckle halo. This is 
useful as a description of the mean starlight intensity, and the derivations are clear and familiar since the product of the 
field and its conjugate lead to phase differences and hence phase structure functions that are insensitive to large-scale 
phase wander. The same is not true of the field, which has an overall phase that possibly wanders over many 2n cycles. 
Even a highly-corrected post- AO field can wander in phase since the wavefront sensor and any corrected images are 
blind to the phase. If we somehow measured the complex halo field directly while this large uncontrolled piston phase 
wander was occurring, the mean halo field would tend to zero. We get around this problem by using a phase reference 
characteristic of the field over the pupil as a whole, <f> re f(t). The phase- referenced halo field is then 



-#„/(t) 



In our interferometric analysis, we compute the wavefront from the WFS measurements and then use an optical model 
and Fourier optics to compute an estimate of the halo field. We tip-tilt stabilize the PSF images and remove the mean 
slopes from the WFS measurements, so the phase of the PSF core is equivalent the halo field on the optical axis. (Note 
that if we were modeling a coronagraph, the corresponding phase reference would be defined by the PSF core before 
encountering the focal plane mask.) The reference phase is therefore 



»ref 



(t) = arg I j e^ x <*)ll(x)d 2 x j (A18) 



which fluctuates somewhat depending on the residual phase pattern. We can now write the mean phase-referenced 
pupil field 

T(x) = / e *M*.*)-*r./(*))\ T (x). (A19) 

We hav e already made use of the average of a phasor with a GRV phase when we found the phase structure functi on in 
Eq. IA51 Since D/£ao 3> lj many statistically- independent regions of the pupil contribute to the integral in Eq. IA181 
giving a result that is at least approximately a GRV. Assuming further that the AO correction is at least fairly good, 
we can assume that a v < n/2 or better. In very high Strehl, the integral is many sigmas from the origin and the 
complex argument approaches a projection, giving a Gaussian distribution for the reference phase. Its scatter will also 
presumably be much less than a v , by a factor of ~ £ao/D based on the number of independent contributing patches. 
Therefore, we can argue that <f> re f(t) is a GRV and uncorrelated with any given <^(x,t) by at least the area factor of 
{£ao/D) 2 . With these assumptions, the mean pupil field becomes 

T(x) = e-^ 2 ) /2 e-^f) /2 T {x) = e-(^/) /2 J^ 1 / 2 T (x). (A20) 

If the Strehl ratio is high, as it is when working with corrected MMT images in the mid-infrared, we can neglect the 
reference phase variance, giving us the mean halo formula we use in our analysis 

T(x) = ^/a-ToM. ( A21 ) 
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Carrying out the same phase referencing and averaging on Eq. ED with high-Strehl average pupil field, Eq. IA211 we 
can inverse Fourier transform to find the mean pupil field from the mean (or static) phase-referenced halo field 

T oWn(x) = J j^pqi ■ (A22) 

Therefore, once we have an estimate of the static halo, we can use Eq. IA22l to estimate the actual pupil field, which 
includes the NCP aberrations. Since the Strehl ratio is real, the NCP phase aberrations are independent of the Strehl 
ratio and are given by 

^cp(x) = arg U e- iK *V(K)d 2 K \ . (A23) 

The pupil field is the mean plus the zero-mean residuals, which contribute to the focal plane speckles as described 
above. The speckle field at any given instant is the sum of the complex contributions from each statistically independent 
patch, rotated in the complex plane by the local Fourier kernel. If D/£ao 3> 1, there are many independent patches, 
leading to zero-mean complex Gaussian speckle fields, even though the pupil field is far from Gaussian. 

B. GAUSSIAN COMPLEX RANDOM VARIABLES 

We use some less familiar properties of complex Gaussian random variables in this paper and to assist the reader, 
we will provide some derivations. We consider an ensemble of complex values ip n G C with a Gaussian, or "normal" 
distribution. The Gaussian distribution has only two independent moments: the mean /i = (ip) and the standard 

/ 2\ 1/2 

deviation a = ( [i/j — (ip)\ ) . For simplicity, we will limit the discussion here to isotropic distributions (i.e. where 

the standard deviation depends on direction in the complex plane), but the properties we require generalize to non- 
isotropic distributions as well. The isotropic probability distribution is 



e 



-\>l)-H\ % /2o 



oo poo 



Pr(V) = ^ ■ ( B1 ) 

As usual, if i/> — x + *C = P e ^ + P > (X; Pi € K, the probability distribution is normalized 

Pr(V)d X dC = 1 (B2) 

and the ensemble average of a quantity is computed by integrating the quantity weighted by the probability distribution 

poo POO 

(f(4>))= / / /(VOPr(V)dxdC (B3) 



— oo J — oo 



— oo J — oo 



This gives the mean as 

fi= I I ^Pr(^)d X dC (B4) 



oo poo 



oo J — oo 



and the variance 

a 2 = I i ^-Ml'PrMdxdC. (Br,, 



OO fOQ 



— oo J — oo 



The lesser known results that we require are Y = (ip 2 ) (as opposed to the more important X = (ipip*)) and the third 
moment Q = (ip 2 ^)*). For our purposes, -0 is zero-mean, which somewhat simplifies the derivations here. We wish to 
prove that the ensemble average of both Y and Q are zero, and that an A^-sample estimate of each has a scatter that 
is proportio nal to l/\/~N. 

Using eq. IB3I we can write 



oo />oo 



Y = (^ 2 ) = / ^ PrWdxdC (B6) 



— oo J — oo 



v r -p 2 /2a 2 

p 2 e 2Me -—^ P &Mp (B7) 

! d " L" ^' 1 "- (B8) 



Therefore 
Similarly, 
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OO POO 



-oo J — OO 
OO />7T 



i) 2 %l>* Pr(V')dxdC 

« e -id , 



27TCT 2 



-dp 



27TCT 2 



=0 



Q = 0. 



25 
(B9) 

(BIO) 
(Bll) 
(B12) 

(B13) 



Therefore, for zero-mean isotropic Gaussian distributions, both moments are indeed zero. 

We also wish to know how the terms converge to zero in the discrete case with N measurements. An experiment 
producing N samples of ip allows us to make an estimate of the Y moment by averaging the squares of the individual 
measurements 



(B14) 



We can calculate the scatter in this measurement by considering an ensemble of the above iV -mea surement experiments, 
and computing the variance of Yat. The ensemble average of Yjv is obviously zero from Eq. IB9I However, the variance 
of Y/v is not zero and is given by 

(B15) 
(B16) 



t 4 = (\Y N -<y N )\*) = (\Y N \' 

2, 



1 



N 



E^ 2 



N N 



If all of the samples tp n and ip m in a given experiment are statistically independent and zero mean, then for n ^ m, the 
ensemble average of the "off-diagonal" terms are zero since (V'nV'm ) — yPn) yl'm) an< ^ each of the ensemble average 
Y moment factors are from Eq. IB9j) . The "diagonal" terms where n = m are the same sample and do not average to 
zero. We continue 



\n— 1 
1 N 



4(1*1' 

where we have assumed that all of the measurements have the same statistical moments. Therefore 

1/2 



(hAI 4 ) /Vn. 



Similarly, an iV-term estimate of Q is 



The variance of Qn is 



1 N 

q n = (^ 2 r) N = ^ E v£< 



<^q = \\Qn - (Qn) 

N 

E^; 



10 



N 



1 



n=l 
N N 



(B17) 

(B18) 
(B19) 

(B20) 

(B21) 

(B22) 
(B23) 

(B24) 
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Once again, so long as the measurements for n ^ m are statistically independent, so the off-diagonal terms vanish due 
to Eq. IB13I The diagonal terms remain giving us 

i n 

n=l 

Therefore 

<Jq = (M 6 ) 1/2 /^V- (B25) 

The same convergence result holds for any "unbalanced" moment where the number of straight and conjugate factors 
are unequal. 
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